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In a statistical cluster or loop model such as percolation, or more generally the Potts models or 
O(n) models, a pinch point is a single bulk point where several distinct clusters or loops touch. In 
a polygon V harboring such a model in its interior and with 2N sides exhibiting free/fixed side- 
alternating boundary conditions, boundary clusters anchor to the fixed sides of V ■ At the critical 
point and in the continuum limit, the density (i.e., frequency of occurrence) of pinch-points between 
s distinct boundary clusters at a bulk point w £ V is proportional to 

(V>l(wi)^l(W2) • ■ ■1pl{W2N-l)lpl(W2N)^s(w, V)))-p. 

The Wi are the vertices of V, ipi is a conformal field theory (CFT) corner one-leg operator, and ty 3 is 
a CFT bulk 2s-leg operator. In this article, we use the Coulomb gas formalism to construct explicit 
contour integral formulas for these correlation functions and thereby calculate the density of various 
pinch-point configurations at arbitrary points in the rectangle, in the hexagon, and for the case 
s = N, in the 2iV-sided polygon at the system's critical point. Explicit formulas for these results 
are given in terms of algebraic functions or integrals of algebraic functions, particularly Lauricella 
functions. In critical percolation, the result for s = N — 2 gives the density of red bonds between 
boundary clusters (in the continuum limit) inside a rectangle. We compare our results with high- 
precision simulations of critical percolation and Ising FK clusters in a rectangle of aspect ratio two 
and in a regular hexagon and find very good agreement. 

Keywords: pinch point, red bond, conformal field theory, Coulomb gas 



I. INTRODUCTION 

We consider critical bond percolation on a very fine square lattice inside a rectangle 1Z with wired left and right 
sides. Of intrinsic interest to the system are bonds whose activation or deactivation will respectively join or disconnect 
the percolation boundary cluster anchored to the left side of 1Z from that anchored to the right side. Such a bond 
that connects them is an example of a red bond [1]. Red bonds inherit their name from the following scenario. If we 
suppose that only activated bonds conduct electricity and that the wired left and right sides of 1Z are attached to 
opposite leads of a battery, then an activated red bond carries the total current and is hottest, and its deactivation 
stops the flow of current. Red bonds carry similar significance in other physical scenarios modeled by percolation. 
Many of their properties have been studied before, first in context with cluster ramification [2]. The average number 
of red bonds weighted by cluster size was measured in [1] , and the fractal dimension of the set of red bonds is predicted 
in [3, 4], and measured in [5]. Further fragmentation properties of percolation clusters are considered in [6]. In this 
article, we calculate the density (i.e., frequency of occurrence) of red bonds at a given bulk (i.e., interior) point weR 
and some generalizations which we now explore. 

In percolation, red bonds are marked by pinch points, or bulk points where distinct percolation clusters touch. 
We consider the two boundary arcs (i.e., perimeters of the boundary clusters) of the percolating system in 1Z. At 
the center w of a red bond, the two boundary arcs pass very close to each other, separated there by only the red 
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FIG. 1: A percolation configuration with the red bonds, or two-pinch points (colored red) and the boundary arcs connecting 
the vertices of the rectangle pairwise (colored green and blue). The left illustration is a sample in the discrete setting while the 
right figure shows only the (filled) boundary clusters of a sample in the continuum limit. 

bond (figure 1). In the continuum limit, four distinct boundary arcs appear to emanate from w, each ending at a 
different vertex of 1Z. In reality, these four curves are not distinct but join pairwise at (or very close to) w to form 
two boundary arcs. Each boundary cluster is pinched into a narrow channel between an adjacent pair of boundary 
arcs, and they touch each other at (or pass very close to) w where the tips of these channels meet (or almost meet). 
Thus, we call w a pinch point [7] (figure 1). The detail of whether or not the red bond at w is activated is lost in the 
continuum limit where, formally speaking, bonds do not exist but their clusters do. But the location of the red bond 
remains. It is marked by the pinch point at w. Thus, the continuum limit of the density of red bonds in 1Z equals the 
density of pinch points between the left and right boundary clusters in 7Z. 

The connection between red bonds and pinch points generalizes the problem of computing the red bond density 
in TZ to computing the pinch-point density involving s boundary clusters in a 27V-sided polygon V . In particular, 
we suppose that V harbors critical percolation in its interior and exhibits a free/fixed side-alternating boundary 
condition (ffbc). That is, the boundary condition (be) of the sides of V alternate from fixed, or "wired," (i.e., all 
bonds activated) to free (i.e., no conditioning imposed on the activation of the bonds). In this article, we label an 
ffbc event with the symbol In such an event, the ffbc conditions a boundary cluster to anchor to each wired side. 
Now, we define an s-pinch point to be a bulk point w £ V where s distinct boundary clusters touch (or pass very 
close). In the continuum limit, 2s boundary arcs emanate from w, each ending at a different vertex of V (figure 1). 
Clearly, we must have 1 < s < N since at most N distinct boundary clusters can anchor to the fixed sides of V. When 
s = 1 we define a one-pinch point to be a bulk point touched (or approached) by just one of the boundary arcs. As 
the continuum limit is approached, the density of pinch-point events decays as a power law of the shrinking lattice 
spacing (section IIIB). This power, with other scaling exponents, is determined in [3, 7). (For 2A^-sided polygons with 
N > 2, the density of red bonds is still dominated by pinch points involving two clusters in the large system limit 
since s-pinch points with s > 2 occur much less often, as discussed in section IV.) 

We obtain another generalization by considering the statistics of the boundary arcs, which fluctuate in V with the 
law of multiple- SLE K [8]. The case of percolation entails K = 6, but we may consider other k £ (0,8) as well. In 
these terms, an s-pinch point is a bulk point w £ V where s distinct multiple-SLE K curves touch (or pass very near 
each other). In particular, a one-pinch point is a bulk point on (or very near) one of these curves (figure 2), and the 
problem of calculating its density generalizes the same problem for when there is one SLE K curve. The latter was 
originally solved in [9]. In our situation with multiple boundary arcs, the regions that a boundary arc can explore in 
V are influenced by the presence of the other boundary arcs, so a one-pinch point can be interpreted as measuring 
the repulsion between the various boundary arcs. In the case of percolation, this "repulsion" is not felt until the 
boundary arcs actually collide due to the locality property of SLE K with n = 6 [10, 11]. 

The range k £ (0, 8) describes boundary arcs in many interesting critical lattice models, including those of the 
Q-state Potts model for Q < 4. As in percolation, an s-pinch point is still a bulk point where s distinct boundary 
clusters touch in the Q-state Potts model inside V with an ffbc. However, now there are two different types of clusters 
to consider: FK clusters and spin clusters. Boundary arcs of the former type are multiple-SLE K curves with speed k 
in the dense phase (i.e., k £ (4, 8)) and related to Q through [12] 

Q = 4cos 2 (47t/k), kg (4,8), (1) 

and boundary arcs of the latter type are multiplc-SLE K curves with "dual" speed k = 16/ K in the dilute phase (i.e., 
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FIG. 2: An illustration of one-pinch-point events on the perimeters of the boundary clusters (orange and purple) in the discrete 
(left) and continuum (right) settings. 

k £ (0,4]) and k related to Q through (1). Scaling exponents and fractal dimensions associated with pinch points are 
found in [4] . The generalization of red bonds from percolation to other models is also considered in [5] . 

In this article, we calculate various continuum limit pinch-point densities in the rectangle 1Z and in the hexagon 
T-L (and for s = N in any 2A^-sided polygon) conditioned on a specified ffbc event ? and for arbitrary n £ (0,8), 
but before we begin, we refine our definition of a pinch point. We suppose that a multiple-SLE K process evolves 2N 
distinct boundary arcs anchored to the vertices of V until they join pairwise in the long-time limit to form N distinct 
boundary arcs in one of C/v possible connectivities. Here, Cjv is the N-th Catalan number, given by 



(27V)! 



N\(N + 1)\ 



(2) 



We let A label a pinch-point event, that is, an event containing all boundary arc configuration samples in which s 
distinct boundary arcs, each with both endpoints among 2s specified vertices of V, pass within a small distance 5 
from a specified bulk point w £ V and the other boundary arcs join the remaining vertices of V in some particular 
connectivity. Then for a specified ffbc event <;, the type-A s -pinch-point density pT^\{w) is the probability of the 
pinch-point event A conditioned on the ffbc event c, and it equals the ratio of the (continuum limit) partition function 
^(Ak) summing exclusively over samples in A fl <? divided by the (continuum limit) partition function Zf summing 
exclusively over samples in 

The purpose of this article is to study the asymptotic behavior of the type-A pinch-point density as S — > 0. The 
asymptotic behavior of the partition functions ^T\u\ arL d Z^ are supposed to be 
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(3) 



where Zj is the free partition function (summing over all samples in the system configuration space), where the 
functions T % and are universal partition functions, where 9\ is the boundary one-leg weight associated with 
the free/fixed boundary condition change (bcc) at each vertex of V , and where S is the bulk Is-leg weight associated 
with the bulk s-pinch-point event. Also, c\ is a non-universal scaling coefficient associated with each bcc, C s is a 
non- universal scaling coefficient associated with the s-pinch-point event (and is not the s-th Catalan number), and 
the i-th bcc occurs within distance Si from the i-th vertex Wi of V . Then the density behaves as 



= Z, 



(A| 



~ C s 2 5 2e =Tf A| 

<5->0 



S )/ T 



(4) 



Thus, determining the behavior of pT\i^ 



to within a constant amounts to determining the universal partition functions 



Tf A|f) and T*. 

The organization of this article is as follows. In section II, we identify the universal partition function T j^, , with a 
bulk-boundary CFT correlation function of certain primary operators, and we find an explicit formula for it using the 
Coulomb gas formalism. Also in this section, we calculate the A^-pinch-point weight (defined below) of a 2Af-sided 
polygon, and we find that it is completely algebraic. In section III, we compute various s-pinch-point densities in the 
rectangle (N = 2) and in the hexagon (N = 3). We find that the formulas for the (N — l)-pinch-point densities are 
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given by algebraic factors times Lauricella functions of cross-ratios of (the half-plane conformal images of) the bulk 
point w and the vertices Wi of V. In section IV, we compare some of our predictions with high-precision simulations 
of percolation and Ising FK clusters inside a rectangle and a regular hexagon and find very good agreement. 



II. CONFORMAL FIELD THEORY DESCRIPTION 

In the continuum limit, l"^^ equals a correlation function of appropriate primary operators belonging to a con- 
formal field theory (CFT) [13] of central charge [14] 

c = c(«) := (6- k)(3k-8)/2k. (5) 

These primary operators are chosen as follows. To alternate the be on the sides of V from free to fixed to free, etc., 
we insert a corner one-leg operator ip1{wi) at each vertex Wi of V into the correlation function. Corner operators are 
defined in [15, 16], and they are used in section IIIC. The collection of 2N corner one-leg operators introduces N 
non-crossing boundary arcs that connect the vertices w%, . . . , w 2 n pairwise in one of Cn possible connectivities. Now, 
to generate an s-pinch point at w € V, we require s of these arcs to touch at (or come very close to) this point. One 
may view this as the event in which 2s distinct boundary arcs emanate from w, which is conditioned by the insertion 
of a spinless bulk 2s-leg operator ^ s (w,w) into the correlation function [7, 17] (figure 3). Hence, T(ai c ) is given by 
the (2N + 2)-point function 

t Fa| ? ) = (^i(wi)ipt{w 2 ) ■ ■ .iPi(w 2 n-i)^i(w2n)^s(w,w)) v . (6) 

The standard approach to studying this correlation function is to conformally map the interior of V onto the upper 
half-plane EL The half-plane version of this correlation function is 

T (A|?) = {%i)ipi (x 2 ) ■■ -ipi (x2N-i)^i(x2n)^s(z, z))h (7) 

= (lpl(xi)lpi(x 2 ) . . .1pi(x 2 N~l)lpl(x2N)^s(z)^s{z))c, (8) 

where we have used Cardy's method of images [15] to rewrite the half-plane correlation function on the right side of 
(7) as the whole-plane correlation function (8). Here, ^ s ( z ) has holomorphic weight Q s and antiholomorphic weight 
zero. First, we will focus on calculating T(a|c)j and later in section IIIC, we will transform T(A| f ) to "f^,^. 

In the multiple-SLE K picture, the bulk 2s-leg operator conditions a specified collection of 2s of the 2N available 
multiplc-SLE K curves to grow from their respective origin points at the vertices of V towards the common bulk 
point w G V until they join pairwise very near w in any one of C s possible connectivities. Here, C s is the s-th 
Catalan number (2). The s-pinch-point event A, defined above, contains all samples that exhibit any one of these C s 
connectivities near the pinch-point. Each connectivity is equally likely to occur. 

Moreover, we may view each pinch-point sample in A as a collection of loops that surround the perimeters of the 
bulk and boundary clusters in V . These clusters are, for example, FK clusters or spin clusters in a Potts model, and 
the loops for the former case are shown in figure 4. This picture is consistent with the continuum limit of the O(n) 
model if we set the loop fugacity n equal to [10, 18] 

n = n(n) := -2 cos(47t/k). (9) 




FIG. 3: An s-pinch-point event is induced by the insertion of a bulk 2s-leg operator. 



FIG. 4: The unique loop configuration corresponding to a particular bond configuration in 1Z with the left/right sides wired. 
Boundary loops (surrounding boundary clusters) are red while bulk loops (surrounding bulk clusters) are any color but red. 



In this interpretation, the endpoints of the boundary arcs are joined pairwise via N exterior arcs to form between one 
and N boundary loops (red loops in figure 4) that dodge in an out of V [16, 19]. These exterior arcs live outside V 
and connect its vertices pairwise, and their connectivity is determined by the choice of ffbc. All samples in the same 
ffbc event ? whose boundary arcs join in the i-th connectivity inside V have the same number pi of boundary loops. 
Because the s-pinch-point event A sums over all C s possible connectivities of the boundary arcs that are conditioned 
to approach the s-pinch-point, we may factor out the fugacity factors associated with the boundary loops to write 



"F(A| f ) m t ne form 



T (Ak) = (n Pl +...+ 



)n A . 



(10) 



The factor ILa is called the (half-plane) type- A pinch-point weight, and it bears the same partition function interpre- 
tation as T(a| ? ), but with the boundary loops having fugacity one. We will elaborate on the relation between an ffbc 
event and the number of boundary loops in each of its samples further in section III B. 

The bulk 2s-leg and boundary s-leg operators ^ s and tj) s respectively are primary operators of a boundary CFT in 
the upper half-plane. The highest-weight vector of their Verma modules respectively belongs to the (0, s) and (1, s + 1) 
(resp. (s,0) and (s + 1, 1)) positions of the Kac table in the dense phase (resp. dilute phase) of SLE K [17]. The Kac 
weights associated with these positions are given by [13] 



h r .s ( K ) 



1 - c(k) 
96~ 



r + s + (r — s)^ 



'25 - c(k) 
1 - c(k) 



1 } (Kr- 4s) 2 - (k - 4) 2 
16k I (ks ~ 4r) 2 - (k - 4) 2 



k > 4 
k < 4 



(11) 



Thus, the conformal weights O s and 9 S of the boundary one-leg operator and the bulk 2s-leg operator in either phase 
are respectively 



16s 2 - (k - 4) 2 
16k 



s(2s 
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2 k 



(12) 



CFT translates the reducibility of the Verma module associated with each boundary one- leg operator into the following 
semi-elliptic system of 2A^ PDEs that govern the correlation function (7), or equivalently, the pinch-point weight Ha: 
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d, 



i X 3 



d z 

Z Xj, 



(z - Xi) 2 



Z - Xi 



e, 

(z - Xi) 2 



n A = o, i e {1,...,2A}. (13) 



The domain of 14a is such that xi < Xj when i < j and z and z are in the upper and lower half-planes H and H* 
respectively. We treat z and z as independent holomorphic and antiholomorphic variables until the very end of our 
calculations, where we set z = z* . (Throughout this article, "z*" denotes the complex conjugate of z.) In addition, 
the three conformal Ward identities ensure that IIa is conformally covariant such that each boundary point x^ has 
scaling weight 6\ and the bulk point z and its image point z have holomorphic weight 9 S : 



2N 



dz + dz + ^2 d Xi 



n A = o, 



(14) 
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zd z + zd- z + 2Q S + yVsiSc, + 9 X ) 



i=l 



n A = o, 



z 2 d z + z 2 d- z 



2N 



29 s (z + z) + }_^{xid Xi + 29-LXi) 



n A = o. 



The Ward identities restrict IIa to a conformally covariant ansatz which may be chosen to be 

N 

n A (a;i, . . . , x 2N ; z, z) = \z- z\~ 20s Y\_( X2i 



X2i-l) Wl G(?72, ... ,V2N-2; M> v), 



(15) 
(16) 

(17) 



where {772, . . . , f?2JV— 2> ^} is a maximal set of independent cross-ratios that can be formed from the points 
X\, . . . , X2N, z, and z, and where G is an unspecified function differentiable in each independent variable. We choose 



(xj - xi){x 2 n - acgjy-i) 

{x 2 N-l - X 1 )(x 2 N - %iY 



(z - Xi)(x 2 N ~ X 2 N-l) 
{X 2 N-1 ~ X 1 )(X 2 N - Z) ' 



(z - X\){x 2 N " ■EgjV-l) 

(a;2jv-i - ^i)(^2jv - 2)' 



(18) 



This ansatz reduces the number of variables in the problem from 2N + 2 to 2N — 1. A standard approach that takes 



advantage of this reduction is to transform (13) into a system of PDEs governing x 2A }IIa and take the limit 

{Xx, X 2 , ■ ■ ■ X 2 N-2,X2N-l,X2N, Z, z} {0, T) 2 , ■ ■ ■ , V2N-2, 1, 00, /It, v} . 



(19) 



This gives a system of 2N PDEs governing the unknown function G from which we can glean information, ideally 
exact solutions. Because we mainly consider the cases N = 1, 2 and 3 in this article, we use the following notation 
throughout: 



V '■= V2, 



V3, & ■= m- 



(20) 



We can explicitly solve (13-16) in the case s — N — 1 (i.e., the two-gon). In this one-pinch-point event, a boundary 
arc 7 connecting x\ with X2 passes some very small distance e from the specified bulk point z € H. We denote the 
corresponding pinch-point weight by II12. Substituting the ansatz 



U l2 {x l ,x 2 - 1 z,z) = \z-z\- 2ei {x 2 -x 1 )-'' ai G{v) 



(xi - z)(x 2 - z) 
(xi - z)(x 2 - z) 



(21) 



(slightly modified from (17)) into (13) yields a second-order, linear, homogeneous differential equation in G. The 
general solution is 



F(xi,x 2 ;z,z) = 



(x 2 - xx) 



-201+02 I 



-2ei+0 2 



\Xi 



\X 2 



;|8/k-1 



ci + c 2 (3 



1 



(xi - z)(x 2 - z) 
{xi - z)(x 2 - z) 



(22) 



where /3 is the incomplete beta function, c\ and c 2 are arbitrary real constants, and the weights #i,#2, and 0i are 
given in (12). We argue that C2 = in our application by sending z — > x € M \ {xi,X2}- Because the boundary arc 
7 is conditioned to touch z, 7 will touch the real axis at x in this limit, and two boundary arcs will emanate from 
x. Thus, the bulk operator ^i(z) must fuse with its image ^i(z) to create a boundary two-leg operator i/j 2 (x) to 
leading order. Or instead we send X2 — > xi. Then in this limit, the two endpoints of 7 touch at x\, and the boundary 
operators ipi(xx) and il)\{x 2 ) fuse to create ^{xi) to leading order as well (figure 5). In both cases, v — > 1. Because 



P(a,b\v) 



if 6 < 0, 



(23) 




FIG. 5: The one-pinch-point configuration in the upper half-plane. We note that the limits X2 
each generate a boundary two-leg operator to leading order. 



a?i and z — > x G R \ {xi, X2} 
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and because b = —62 < for k < 8 (12), we see that 

C2 / (X2 — X\){z — Z) 



{x 2 -x 1 y 20 ^\z-z\- 2Q ^ 

F( Xl ,x 2 ;z,z) ^ \ Xi _ z \ S / K -i\ X2 _ z \8/ K -i 



Cl 



2 Vfcl - ^)(^2 - ^) 



(24) 



To ensure that the bulk-image or boundary-boundary fusion has the two-leg channel at leading order, the second term 
in the brackets must be absent. Thus, c 2 — 0, and we find the one-pinch-point weight for an SLE K connecting the 
real points x% and x 2 : 

(x 2 -x 1 f/-\z^z\^l^ 
^{xux 2]z ,z) = 1 — w -^—^. (25) 

If we put X\ = and x 2 — 00 as in the usual setup for SLE K , then we have 



- 2ei P{S(e,z)n 7 ^0}, «e(0,8), (27) 



lim ^2 til U 12 (0,x 2 ;z,z) = \z - z\ K l % ~ x sinarg(z) 8/K - 1 (26) 

X2— ^OO 



where P{S(e, z) n 7 7^ 0} is the probability that 7 intersects a ball Z?(e, z) centered at z <E H and of small radius e. 
Equation (27) is rigorously proven in [9]. This rigorous result is supposed by physicists to be stronger. Namely, it is 
expected to be 

P{£(e, z) H 7 ^ 0} ~ Ce 2ei |z-z| K / 8 - 1 sinarg(z) 8 / K - 1 , k g (0, 8) (28) 

e— >0 

for some constant C. This is equivalent to the prediction (4) when N = 1. We note that 0i > for k < 8, so this 
probability goes to zero as e — >■ as it must. Below, we will compute this pinch-point weight with another method. 

As we observed in this example, it appears to be generally true that the set of pinch-point densities span a proper 
subspace of the solution space of the system (13-16). This follows from the result (A23) in the appendix A. 

The system (13-16) is very difficult to solve directly when N > 1, but fortunately the Coulomb gas formalism [20] 
provides a tractable approach to constructing explicit solutions. To this end, we write a chiral operator representation 
for (8). That is, we represent a primary field of holomorphic weight h and antiholomorphic weight zero by a chiral 
operator V a (z) of charge a. A primary field, this chiral operator is defined to be the normal ordering of a exp[i\/2a.(p(z)] 
with <p(z) the holomorphic part of a massless free boson, and its holomorphic weight is h — a(a — 2«o)- Here, 2ckq 
is the background charge, and it equals a + + a_ with the screening charges a± given in (29) below. Only chiral 
operators of charge a or 2«o ~ ct have equal holomorphic weights, so we call these two charges conjugates. The two 
conjugate charges af s associated with the Kac weight h r s are 

+ l±r lis , ((Jk/2)* k>4 

^ = — a+ + — a ~> a± = ± ((^/2F «<4' (29) 

and we let Vjr a (z) be a chiral operator of charge a^ s . Adopting dense phase (i.e., n > 4) notation conventions, we 
represent ^i(xj) by the chiral operator V^~ 2 {xi), and we represent ^ s (z,z) by the vertex operator V^ s {z)V^ s {z). The 
correlation function (8) now has total charge 

2Na i.2 + 2a ls = 2a + (s- N)a-. (30) 

We wish for this total charge to equal 2a in order to satisfy the neutrality condition. This is necessary in order 
for the correlation function to satisfy the conformal Ward identities (14-16). We see that the neutrality condition is 
presently satisfied only when N = s. 

We momentarily restrict our attention to the case s = N, where the correlation function is neutral. Here, we find 
an explicit, algebraic formula for the upper half-plane A-pinch-point weight in a 2A-sided polygon: 

njv.pinch P oint(^i, • ■ • ,x 2N ; z, z) = (V{^ 2 (xi) . . . V^ 2 {x 2N )V^ s (z)V£ s (z)V-(ui) . . . V-(u N - s )) 

2N 2N 



\ Z _- Z \(AN+A- 



K) /8K Il^i - X ^ /K II l z Xi| 1 - 4( * r+1)/K . (31) 



i<j 



The bulk point z is connected to all boundary points x\ hosting the becs via the N boundary arcs that touch at z. 
We note that (31) is identical to (25) when N = 1, as it must be. This result was also computed in [21, 22] by using 
other methods. 
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= 4e 7ri(/3l -' 92) sm7r/3 1 sin7r# 2 r"~*~^ » 

%x x 2 

FIG. 6: The Pochhammer contour entwining the points xi and X2- If the numbers /3i and 02, where e 27rl/31 and e 2 '"' 32 are the 
monodromy factors associated with x\ and respectively, are greater than negative one, then a Pochhammer contour may be 
replaced with the simple contour shown on the right. 

Next, we consider the cases with s < N. In order for the total charge (30) of the correlation function to equal 2a>o, 
we must insert N — s Q- screening operators, leading to the following modified (2N + 2)-point function: 



(VT 2 ( Xl ) . . . V^ 2 {x 2N )V+ s (z)V+ s {z)V-{u x ) . . . V-(uff- a )) d Ul ... du N ^ s . (32) 

ri Jr N - e 

After including a useful prefactor discussed below, we find that (32) is given by 



Tf .Mrp-iW A fTT(.-.-H (j\ 

Jj- i 4ex P 7rz(/3 1 „ l -/3 2m )sin^/3 lm sin 7 r/3 2m r(l-4//c) 2 M 11 J l) Mil 



z — x. 



I (re— 4s— 4) / re 



i<j ) \i=l 

J.\ \ •• \ /N-s 

-•■ I ■ - zl^ -48-4 ^ 78 " <Z> ... <P dux... dun-* III II (.<•;. - a, \ '"II I I ! //„ - // 



I ... I dux... du N - s I f[ Yl (x k - Ul )- 4 ^ J ( H (i 

JTi Jr N ^ s \k=l 1=1 J \p<q 

x ( n (z - u m )^ +i -^^(z - u m )^+ 4 -^A . (33) 

\m=X J 

and e 27Tl ^ 2m are the monodromy factors (relative to u m ) of the two branch points entwined by the m-th 
contour, n(n) is the O(n) model loop fugacity (9), and we choose the branch of the logarithm with argz € [— 7r,7r) 
for all complex z so that each branch cut parallels the real axis. Every pinch-point weight Ha will be some linear 
combination of functions of the form (33), with each term using a different set of contours {T m }. A proof that (33) 
solves the system (13-16) is given in section A. 

What remains is to determine a collection of closed, non-intersecting integration contours {r m }^Z^ appropriate for 
a particular type-A pinch-point event. The simplest closed contour along which an integration is nonzero is a closed 
Pochhammer contour entwining only a pair among the branch points xx, ■ ■ ■ , x 2 n, z, and z of the integrand, as shown 
figure 6. Throughout this article, we take each T m to be such a contour. 

Now we explain convenience of the prefactor in (33). A first reason involves the limit x i+1 ~^ %i that sends the 
27V-sided polygon to a (2JV — 2)-sided polygon. After multiplying the half-plane weight IIa by (xi+x ~ Xi) 201 and 
taking this limit, this product goes to either zero or an s-pinch-point weight for the half-plane conformal image 
of a (2N — 2)-sided polygon with vertices sent to xx, ■ ■ ■ ,Xi-x,%i+2i ■ ■ ■ ,x 2 n- Meanwhile, one can show that (33) 
goes to either zero or the same expression except with all factors containing a;j,a;j-|_i, and un- s omitted, the un-s 
integration omitted, a factor of /3(— 4/rt, -4/k) -1 = T(2 — 8/k)/T(1 — 4/k) 2 (with /3(a, b) the Euler beta function) 
and [4sin 2 (47r/K)] _1 omitted, and possibly a fugacity factor n(fc) omitted. So to within a factor of n(n), we retain 
the same normalization as that for the case of the 2iV-sided polygon with N M> N — 1. 

A second reason involves the cases k = 4/o for some positive integer a. At these special values, all of the powers 
in the integrand of (33) are integers. As a result, the m-th Pochhammer contour entwining two branch points, now 
with their respective monodromy factors e 2m ^ lm and e 2m ^ 2m equaling one, disintegrates into two pairs of oppositely- 
oriented loops. One pair surrounds the first point, the other pair surrounds the second, and each integration thus 
gives zero. More precisely, one can show that the (N — s)-fold integral in (33) is 0((k — 4/a) N ~ s ) as K — > A/a with 
a € Z + . Meanwhile, the complete prefactor (with each /3 lm and j3 2m equaling — 4//t or 4(s + !)/« — 1) is 



N-s 

n 

m— 1 



/ cos(47r/«0r(2 - 8/k) \ , , tN sU , , . . 

-— a y / > K ' ^ = O ({k - 4/a- (JY ~ a) as K ^ 4/a, a e Z+ . 34 

V 2sin7r/3i m sm7r^ 2m r(l-4/K)V 



Therefore, the product of the prefactor (34) with the integral in (33) is finite and nonzero in the limit. We note that 
n(«)r(2 — 8/k) is also finite when k — > 8/a with a an odd integer greater than one. These two cases cover all of the 
singularities of the gamma functions appearing in (34) for k € (0, 8). 
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A third reason involves simplifying the contours. In some cases, a Pochhammer contour entwining two branch points 
may be replaced by a simple curve that starts and ends at the those points (figure 6) [23] . Each branch point has some 
monodromy e 27 ™^ relative to each integration variable: for each Xi, ft = — 4/k, and for z and z , ft = 4(s + 1)/k — 1. 
The former power is greater than negative one only when k > 4, so only then can we replace a Pochhammer contour 
entwining an Xi by a simple curve. Otherwise, such a replacement yields a divergent integral. The latter power is 
greater than negative one for all k € (0, 8), so any Pochhammer contour entwining z with z can be replaced by a simple 
curve. Replacement of the m-th contour by a simple curve cancels the factor of 4 exp 7ri(/3 lm — /3 2 m) sin 7r/3 lrra sin 7r/3 2m 
in the denominator of (33). Throughout this article, we will explicitly use simple curves in place of Pochhammer 
contours and omit these factors. If one of these replacements creates a divergent integral, then we implicitly revert 
back to using a Pochhammer contour for that integral, and we include the omitted factor. 

III. CALCULATION OF PINCH-POINT DENSITIES 

In this section, we calculate the explicit formula for the type-A s-pinch-point density P(a| ? ) conditioned on the ffbc 
event ? for various s-pinch-point events A for either a rectangle 1Z (N — 2) or for a hexagon T-L (N = 3). We proceed 
in four steps, as enumerated below. 

1. We compute the half-plane pinch-point weight IIa for the pinch-point event A of interest. 

2. We construct from IIa the universal partition function Y(a|^) that sums exclusively over the event An?. 

3. We transform Y(A| S ) into the universal partition function TT^\ with the appropriate polygon V for its domain. 

4. With V = 1Z or H, we divide by the universal partition function that sums exclusively over the ffbc 
event <; to obtain formulas for pinch-point densities in V . 

A. Half-plane pinch-point weights 

1. The case N = 2, s = 2 and N = 3, s = 3 

First, we consider the two-pinch-point weight n 12 34 for N — 2 and the three-pinch-point weight n 12 3456 for N = 3. 
The subscripts on each weight indicate the indices of the points Xi that are connected to the pinch point by a boundary 
arc, and in this case, these are all of the available vertices. The weights are given by (31) with N — 2 and N = 3 
respectively, and both may be expressed in the covariant form (17) which we will find convenient later. We find that 
the half-plane (s — N — 2)-pinch-point weight is 

IW(zi, ...,x 4 ;z,z) = [{x 2 - x x ){x A - x z )] l - & ' K \z - z |«/8-6/«-i 

x rf'^il - ?7) 2 /> - H 24/K ~ 2 M?7 - Ml ~ ")(1 - M)(l - ^)] 1/2 - 6/k , (35) 
and the half-plane (s = N = 3)-pinch-point weight is 

IWseOi, ...,x 6 ;z,z) = \{x 2 - xi)(x A - x 3 ){x 6 - xg)] 1 " 6 ^^ - z|»/8-ie/«-i 

x [ri(a - T)f K - l [Ta(T - V )(a - - r))(l 
x [fj,f(r] -fJ,)(r}- v)(j - fi)(r - v){a - fx)(a 

where the cross-ratios r),T,cr, fi,f are defined in (18, 20). The correct normalization of these pinch-point weights 
depends on bulk-boundary fusion coefficients, but because it is not needed for our purposes, we ignore it in this 
article. 

Before we calculate s-pinch-point weights for s < N, we comment on the normalizations of these weights too. When 
s < N, virtually all samples in the s-pinch-point event A will have at least one interval {xi,Xi+\) with its endpoints 
mutually connected by a boundary arc that does not pass near the pinch-point. The fugacity of this boundary arc is 
one since we are working with pinch-point weights, so when we send 2^+1 — > xi (after multiplying by {xi + \ — Xi) 1 
first), we must recover an s-pinch-point weight (independent of Xi) in a system with becs at the remaining 2(A^ — 1) 
points on the real axis. Continuing this process until no such intervals remain, we eventually reach an s-pinch-point 
weight in a system with becs at the remaining 2s points on the real axis. The weight of this event is given by (31) 
with N — s. Therefore, the s-pinch-point weights with s < N are normalized so that they equal the s-pinch-point 
weight (31) with N s after this sequence of N — s limits is taken. 



r)(l-a)] 2 / K | M -H 48 / K - 3 
^(l- M )(l-^)] 1 /2-8/^ 



(36) 
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X\ X 2 X 3 X4 




FIG. 7: The pinch-point configuration for Il23:4i. The left (resp. right) figure shows that a boundary two-leg (resp. four-leg) 
operator is generated to leading order when z approaches the intervals (3:1, 3:2), (X2, £3), and (3:3,3:4) (resp. the interval (3:4, xi)). 



2. The case N = 2, s = 1 



Next, we consider one-pinch-point events with N — 2 boundary arcs. Here, one boundary arc 71 connects the 
points Xi and Xj with a bulk point z G H while the other boundary arc 72 connects the remaining points Xk and xi. 
The weight TUj-.M of this event is given by (33) with N = 2 and s = 1. 

The formula for Ily.jy has a single contour integral 1^.^, and the contour is chosen so that the chiral operators 
exhibit specific fusion rules that depend on which vertices are connected to z through 71. For example, we consider 
Il23:4i (figure 7). If we let the bulk point z approach a boundary point x in the segment (11,352), then 71 must touch 
(2:1,2:2) at x in this limit, which is equivalent to placing a boundary two-leg operator if>2(x) there. Therefore the bulk 
operator ^fi(z) must fuse with its image ^i(z) to give ip2(x) to leading order. Now, when the two chiral operators 
Vq^(z) and Vq~ 1 (z) fuse, their product is a boundary chiral operator with charge 2a^i = a^ 3 which carries the weight 
02 of a boundary two-leg operator as required (12). The same is true of the intervals (2:2, 2:3) and (£3, X4). 

Next, we let the bulk point z approach a boundary point x in the interval (2:4,2:1). Because 72 joins x± with 2:4, 
topological considerations show that both 71 and 72 must touch (X4, xi) at x in this limit (figure 7). Therefore the 
leading operator of the ensuing bulk-image fusion must be a boundary four-leg operator. Above, we saw that the 
total charge of the bulk-image pair equals that of a boundary two-leg operator. But if we add the screening charge 



then this total charge becomes 2a, 



0,1 



- = a i5: which is that of a chiral operator with the desired boundary 
four-leg weight #4 (12). The screening charge is pulled in with the bulk-image fusion only if r23 : 4i contracts to a point 
in the process. Thus, F 2 3 : 4i must be a simple curve starting at z and ending at z. 

In order for Il23 : 4i to be a continuous function of z and z, each point of r23 : 4i must reside on the same Riemann 
sheet of the integrand, so r23 : 4i can only cross the real axis through a specific segment (2^,2^+1). (Here, 2:5 :~ x%.) 
This segment must be (2:4, x±) in order to ensure that T23 : 4i contracts to a point when we let z and z approach a point 
in (2:4,2:!). This choice creates another desired effect. In the event of a bulk-image fusion over (2:1,2:4) \ {2:2,2:3}, 
r23:4i does not contract to a point, the screening charge is not drawn in, and an undesired boundary four-leg operator 
in (xi, 2:4), which would contradict the assertions of the previous paragraph, is not produced. 
By cyclically permuting the indices, we find four one-pinch-point weights: 



{ni2:34, n23:4i, n34 : i2, n4i : 23j-. 



(37) 



Each weight is given by (33) with N = 2, s = 1, and r^jy a simple curve connecting z with z and crossing the real 
axis only through (xk,xi). In the formula for each weight, we order the differences in the factors of the integrand so 
that the branch cuts do not intersect Tij : ki and the integrand restricted to Tij±i is therefore a continuous function of 
X\ , . . . , 2:4, z, z, and u := u\. 

It is useful to decompose these one-pinch-point weights (37) into a linear combination of the integrals (with 2:5 :— z 
and 2:6 := z) 



du N 



ie{l,...,6}, 



(38) 



in order to explicitly show that these weights are real (or at least share a common phase, as they must be for physical 
reasons) and to express them in terms of Lauricella functions. The operator "A/"" orders the differences in the factors 
of the integrand so that Ii is real. {I\ is a sum of integrations from xq :— xq to 00 and from —00 to x\.) Because 
arg(z) e [— 7r,7r) for z £ C, the integrand has a branch cut that starts at each xj with j < i (resp. j > i) and points 
leftward (resp. rightward) along the real axis. For simplicity, we momentarily suppose that 2:5 and x§ are real as 
we decompose each one-pinch-point weight into a linear combination of the various Ii times algebraic factors. For 
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FIG. 8: The contour used for the one-pinch-point weight Ili2:34. To facilitate calculation, we at times place z and z at adjacent 
locations xs, xe respectively on the real axis as in the right figure. (In each figure of this article, a blue (resp. orange, resp. red) 
circle marks a point of charge aj~ 2 (resp. s , resp. a_) in the dense phase.) 



example, we can use figure 8 to find the decomposition 



n 



12::S-1 



A 



2i sin 



to - * 5 r /8+8/K - 2 Y[( Xj Xl fi- n(^ - x 4 ) i/2 - 4/k (x 6 - x 4 ) i/2 - 4/k . (39) 



The proportionality constant A will be determined below. Now, to show that the one-pinch-point weights have a 
constant phase that may be adjusted to unity, we seek a basis of integrals for the span of {It,... ,Iq} that are real 
when x$ = z and xq — z. Integrating the screening charge along a contour parallel to and immediately above 
(resp. below) the real axis gives the — (resp. +) branch of the linear relation 



h + e ±47 " /K / 2 + e ±Snt/K h + e ±12m/K I 4 + e ±167rz/K I 5 - e ±87ri/K I 6 = 0, 
which allows us to write Jg and Iq in terms of Ii, . . . , I 4 . Then because the integral 



(40) 



Xi := £(-4/k, -4/k)" 1 / duM 



fa _ u f^\x & - ufl K ~ l ]J(u - x z 

3 = 1 



-4/k 



h ~ e M/K I 6 + h (41) 



is real when X5 = z and x§ = z, we find a real basis {Zi,22 := 12,^3 := ^3,^4 := ^4} for the span of We 
anticipate that the coefficients found from decomposing the Hij-.kl over this basis will share a common phase. 

However, it is more useful for our purposes (of calculating one-pinch-point weights for the hexagon later) to compute 
this decomposition via a different approach in which the four integrals Z$ arise naturally as conformal blocks. We 
consider the one-pinch-point weight TLi2(xi,X2',z,z), given by (ipi(xi)\2]' l Pi( x 2)^i(z)[2}'^i(z)) ■ Here, the bracketed 
subscript between a pair of primary fields indicates the unique fusion channel propagating between that pair, so "[s]" 
indicates the s-leg channel when s > and the identity channel when s — 0. To increase N from one to two, we insert 
the charge-neutral collection J^ 4 duVi 2 { x 3)^i2( x i)^-( u ) with x 2 < x$ < X4 into the chiral operator representation 
of this four-point function. We find 



(V'i(2 ; i)[2]'0i(2 ; 2)V'i(2 ; 3)[o]V'i(2 : 4)*i(2)[2]*i(^)) = nJ{x x , . . . x 4 ; z, z) Xi(xi, . ..x 4 ;z,z), 
where the function J is given by (33) with N = 2 and s = 1: 



(42) 



J( Xl ,...,x 4 ;z,z) :=\z-z\ K / 8 +V K - 2 l[(o 



\2/« 



n 



Z Xi 



(43) 



l<3 



The new pair of boundary one-leg operators at 2:3 and x 4 fuse through only the identity channel because the screening 
charge is integrated along a simple curve connecting X3 with x 4 . The original boundary one- leg operators at X\ and 
X2 still fuse through only the two-leg channel. Three of the four one-pinch-point events (ij:kl) are consistent with 
these fusion rules, (41:23), (12:34), and (23:41), so (42) must be a linear combination of the pinch-point weights 
1112:34) Il4i : 23, and Il23:4i. Indeed, this linear combination is (figure 9) 



n 4 l:23 + ^1112:34 + n 2 3 :41 =7lJ(x 1 ,... X 4 ; Z, z)l 4 (x 1 ,. . . X 4 \ Z, z) . 



(44) 



As usual, n is the loop fugacity (9) of the 0(n) model. 

The coefficients of the linear combination on the left side of (44) are found in the following way. First, to find the 
coefficient of n 12:3 4, we send x 4 — > x 3 on both sides of (44). (We always implicitly multiply by (x i+ i — Xt) 2 ® 1 before 
sending Xi+i — > Xi so that the limit exists.) Then Hn^, n23:4i — > while Hi 2: 34 — > H12 (25), and nJ xl 4 -> 71II12. 




FIG. 9: The decomposition of (42) into a linear combination of the weights IIi2:34, Il4i:23, and Il23:4i as given in (44). 



This justifies the coefficient of n that dresses n 12: 34 in (44). Next, to find the coefficient of n 41:23 , we send x$ —> x 2 - 
On the left side of (44), ITi2 : 34, n 2 3 :4 i — > while n 4 i :2 3 — > II14, or really lli 2 with x 2 i-> X4. On the right side, we 
use (40) to write nJ x I4 as a linear combination of nJ x T^nJ x I^,nJ x I 5 , and nJ x I 6 . (These are the four 
nJ x Ii that have either both or neither bounds of integration among {x2,xs}. Again, x^ — z and x§ = z.) All of 
the integrals in this combination except nJ x 7 3 vanish in this limit, while nJ x I 3 goes to rtll 12 . Because nJ x / 3 
carries a coefficient of rT 1 in this linear combination, the right side of (44) becomes IIi 2 with x 2 H> X4. This justifies 
the coefficient of one that dresses H4i :2 3 in (44). The same reasoning gives the coefficient of one for n 2 3 : 4i in (44). 

Cyclically permuting the indices in (44) gives three more equations relating the four one-pinch-point weights (37) 
to the four integrals Zj. Upon inverting these equations to isolate the weights, we find 



n 



ij:kl 



j 



2Xj + (n 2 — 2)li — nli — nlk 



(45) 



n 2 - 4 

6 / K ~ 1 (x4 — X3) 6 / K_1 |z — 2| 1-K / 8 to arrive with a function Gi 



For each index i, we can multiply J x Ti by (x 2 — xi)" /r "~'(X4 — x% 
that depends only on cross-ratios 77, /z, and v, according to (17). After making the replacement (x±, cc 2 , £3, 2:4, z, z) 
(0, 77, 1, 00, /j, 1/), we find 



{ri\n-v\fl K - l {l-r]f/ K 



{nv{q - n)(<q - - /i)(l - i/))V«-i/2 
(HM- ^| 2 ) 4/k_1/2 (1 -r?) 2 / K 

((7,-^(7? -^(1-^(1 - V )yi«-V2 

(r? 2 (l-/x)(l-^)|M-^| 2 ) 4/K - 1/2 



4 4 8 

1 — ; -,1 — , 1 

K K K 



-;2- - 



F 



l-i;*,!-!,!-!;,-! 

UK; Hli hi hi 



1 - 77, 1 - /J, 1 
77 77 s 

/i ' 7/ 



(^(77 - n)(n - 7/)) 4 /«-i/2(i _ v )6/k-i 

(t7|/j-^|) 8 / k - 1 (1 -77) 2 / K 



4 4 
Fd l--;-,l 

K K 



-,1 



-;2 



1-77, 



1 — 77 1 — 77 



FdI 1 



4 4 



1 



1 



(^-/^)^-^)(i-m)(i-^)) 4/k - 1/2 

We have expressed each Gi in terms of the Lauricella function Fjj, defined as [24 

r(o) 



1 -/x' 1 



F D (a,bi 



i ; C I Xi ; . . . , X m ) 



^l-t) - "^!-^) 



-61 



(1 - x m t) 



-b„ 



df. 



(46) 



T(c)T(c-a) J 

by writing the integration variable u of 2^ as the following Mobius transformation of the integration variable t in (46) 

t- 1 



t 



l-(l-ti)t, 



77 1, 
1 

7' 



(47) 
(48) 



These transformations are chosen so that each Fo has m — 3 arguments with the first argument between zero and 
one and the last two arguments being complex conjugates of each other. These choices ensure that each Fjj is real. 
Thus, the half-plane pinch-point weights Hij-.kU expressed in the covariant form (17), are given by 



H ijM {x 1 ,...,x i ;z,z) = [{x2-x 1 )(x i -x 3 )] 1 6/k \z-z\ k/8 1 



2Gj + (n 2 - 2)Gj - nGi - nG k 



(77, /a, v), (49) 
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with r),fj,,v, and n defined in (18, 20, 9) respectively. We note that our normalization in (49) ensures that {x\ — 
Xk) Hij-.kl —> Ilij (25) as xi — > Xk.. Now, if we let k approach A/a with a 6 Z + so that n — > ±2, then the limit 
n.y remains finite. Thus (49) must be finite when n = ±2, although showing this explicitly appears to be difficult. 
Comparing (49) with (39), we find that A — —in/y/A — n 2 . 



3. The case N = 3, s = 2 



Next, we consider two-pinch-point events with N = 3 boundary arcs. Here, two boundary arcs ji and 72, with 
cndpoints respectively at Xi,Xj £ K and Xk,xi G K, touch at a bulk point 2, and the remaining boundary arc 73 has 
cndpoints at x m and x n . We note that this setup restricts the allowed boundary arc connectivities to those in which 
73 does not separate 71 from 72, so x m and x n must be either adjacent or among {x§, x\\. The weight Uijklmn of this 
event is given by (33) with N — 3 and s — 2. By cyclically permuting the indices, we find six such two-pinch-point 
configurations. 

The formula for Hijki-.mn contains a single contour integral Tijki-.mn that is determined via the same reasoning 
that was used for the case N = 2 and s = 1 above. We summarize the argument. The two-pinch-point event is 
conditioned by the insertion of a bulk four-leg operator ^2(2). Topological considerations (as can be understood upon 
examining figure 10) show that fusing this operator with its image across any interval (x a ,Xb) with (a, b) 7^ (m,n) 
(resp. (a, b) — (m,n)) must, to leading order, give rise to a boundary four-leg (resp. six-leg) operator ^4 (resp. ip6)- 



Because 2a„ 



a 



1,0 



is the charge of a chiral operator with the boundary four- leg weight # 4 , this requirement is 



already satisfied when (a, 6) 7^ (m,n). If ^ijki-.mn is a simple curve with endpoints at z and z and crossing the real 
axis only through (x m ,x n ), then the screening charge is drawn into a bulk-image fusion across this interval, shifting 
the product to a chiral operator with net charge 2oq 2 + a - = a t7- This operator has the desired boundary six-leg 
weight 6» 6 (12). 

To express the two-pinch-point weights in terms of Lauricella functions, we write them as linear combinations of 
the six real integrals 



Id :=0(-4/k,-4/k)- 



duN 



(z — u) 



(z-u) 



12/ K- 



1 n>~ x j 

j=i 



.-4/ K 



1 = 1, 



(50) 



(As before, the operator "A/"" orders the differences in the factors of the integrand so that the integrand is real, 
and /Ci is integrated from xq := x§ to 00 and then from —00 to x\.) To proceed, we consider the two-pinch-point 
weight II1234. Inserting the charge-neutral collection duVi 2 { x b)Vi2{ x 5)V-( u ) m t° the chiral representation of its 
six-point function ( - 0i(a ; i)[2]' ! / ; i( a:: 2)'0i( a; 3)[2]' i /'i( a; 4) v f'i( 2: )[4]^ , i(^)) with X4 < x$ < xq, we get the conformal block 



{^lixij^ipi (x 2 )ipi (x 3 )[2]ipi (x 4 )ipi (ac5)[o]V'i( af 6)*i(«)[4]*i ( z )) = nL(xi, . . . , x 6 ; z, z) K&{x\, . . . , x 6 ; z, z), 
where the pre-factor L is given by (33) with N = 3 and s = 2: 

L(x 1: ...,x 6 ;z, z) := \z - s|»/8+«/«-3 f[ {x . _ x ,f/« f[ {z _ ^i-i*/- 

i<j i=l 

After following the reasoning that led to (44), we find (figure 10) 

n 6 123:45 + «n i23 4:56 + n 23 45:61 = nL(Xi, ...,X 6 ; Z,z)K, 6 (x 1 , . . .,X 6 ;Z,z). 



(51) 



(52) 



(53) 



x\ 



+ n 



z O 



■Z'2 



■'"4 



■i-a 



X6 





FIG. 10: The decomposition of (51) into a linear combination of the weights Il6i23:45, Ili234:56, and 112345:61, as given in (53). 
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Another five equations relating the six weights with the six integrals Ki is found by cyclically permuting the indices 
in (53). These equations may be simultaneously solved to give 

'(2 - n 2 ){Ki + JCm) + n{K,j + K,i)-2K, k + (n 3 - 3n)/C n 



n 



ijkl:mn 



= TlL 



(n 2 -4)(n 2 -l) J ' ^ 

(We note the present double-use of n as an index and as the loop fugacity of the O(n) model.) To finish, we seek 
a form for the weights that expresses the /Q in terms of Lauricella functions and exhibits the conformally covariant 
ansatz of (17). To this end, we define the function 

Hi(xi,...,x 6 ;z,z) := [(x 2 - Xi)(x4, - x 3 )(x e - x 5 )] 201 \z - z\ 2e>2 L x ]Q (55) 

= [{xi-xiKxi-xztf' "-^z-z^ 6 ' L 1 xK'i, (56) 

with L and /Cj adjusted to respective quantities L' and K,[ that are finite in the limit x§ — > oo: 

L'(x 1 ,...,x 6 ;z i z) := (x e - x^f'^L, (57) 
lC! i (xi,...,X6;z,z) := (x 6 - x 5 ) 4/K /Ci. (58) 

According to (17), Hi is strictly a function of the cross-ratios 77, r, cr, /i, and v. After making the replacement 
[x\, a?2, £3, Xa, x 5 , xa, z, z) h-» (0, 77, p, a, 1, 00, /j, v), we find 

Hi( V , r, a, fi, v) = Kfa, r, a, /*, u)[r,(a ~ r)] 8 /^ 1 [ra(r rj)(l - - r)(l - cr)] 2 / K 

x l/x - H 24/K - 2 [^(M - - »?)0* - ~ r){t* - a)(u - - l){v l)] 1 /*-*/", (59) 

with each K'^rj, t, a, fi, v) := /C.^(0, 77, r, cr, 1, 00; /J, v) finite and equaling a Lauricella function Fjj times algebraic 
prefactors: 



1 — 77, 1 — T, 1 — (T, 1— (X, 1 — Z/j, 



Jt£fa, r, a, M , v) oc ^^'^{r ^^(a - ^^(1 rj)-^(fi rj) 12 ^ 1 ^ - t?) 12 ^ 1 



37 77 77 77 

T U /i ' ^ 



1 - 



77 cr(r — rj) t — rj fi(r — 77) !/(r — 77) 



r ' r(cr — rj) ' t(1 — 77) ' t(/i — 77) ' t(v — rj) < 
^(77, r, a, M , 1/) oc t" 4 /"(t - Ty) 1 - 8 /^ - r?) 4 ^" 1 ^ - t) 1 ^'^! - r)- 4 / K (/i - r) 12 l^\ V - r) 12 /^ 1 

a - — r i](a — t) [1 — rj){a — t) (/j — rj){o — t) [y — rj){a — r) s 



f d[{X]} (t 77 ' t{o 77) (1-t)(ct-t?)' (^-t)(o--j?)' {v -t)[o - rj))'' 
Kfa, r, a, v) oc (1 r,)- 4 / K (l - r)- 4 / K (l - a) 1 - 8 /^ - M ) 12 /^ 4 (l - t,) 12 /^ 1 

1 — (7 1 — cr 1 — cr 1 — cr 
1 ~ °"' 1 ' 1 ' : i ' 1 

1—77 1 — r 1 — /_i 1 — 7^ 
77, r, a, /J, 7/ 



(60) 
(61) 

(62) 

(63) 

(64) 
(65) 



Again, we have expressed each K,[ in terms of by writing its integration variable u as the following Mobius 
transformation of the integration variable t in (46): 



7 = 1: u — 
i = 3: u = 



t - 1 



T — (t — 77)7; ' 

i = 5: u = 1 — (1 — a)t, 



i = 2: u = rjt, 

rj{a — r)t — r(<r — 77) 
(cr - r)t - (cr - 77) 



i = 4 



7 = 6: u = 



(66) 
(67) 
(68) 



The transformations are chosen so that the first three arguments of each Fd is between zero and one and the last two 
arguments are complex conjugates. These choices ensure that each Fd is real. Each Fd uses the same set of seven 
parameters: 



, ,- 1 4 4 4 4 12 12 8 

{xj}U= !- -;2-- 

■ ' I is* is* is* is* ly* hi* is 1 

rx t% rv r\> r\ 1% r\ 



(69) 
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Combining (54) and (55), we find that each weight is given by the conformally covariant formula 



n 



ijkl:mn (^1 ) 



,x 6 ;z t z) = {(x 2 - Xl )(x 4 - x 3 )(x 6 - xs)} 1 - 6 /^ - z\ K ^ K ~ l 

~n(2 - n?)(Hi + H m ) + n 2 H 3 - 2nH k + n 2 Hi + n 2 (n 2 - 3)H n 
{n 2 -A){n 2 - 1) 



(r],T,a,^,v), (70) 



with each Hi explicitly given among (59-65), r], r, a, fj., and v given by (18, 20), and n given by (9). We note that our 
normalization in (54) ensures that (x n — x m ) 2ei Hijki:mn ~> ^ijki (35) as x n — > x m . If we let K approach the zeros 
4/a and 12/(3a ± 1), with a 6 Z + , of the denominator of (70) in this relation, then the limit Hijki remains finite. 
Therefore, (70) must be finite when n = ±2, ±1, although this seems to be very difficult to prove directly. 

4- The case N = 3, s = 1 

Last, we consider one-pinch-point events with TV = 3 boundary arcs. Here, a boundary arc 71 connects Xi,Xj and 
z, another 72 connects x k and xi, and the last 73 connects x m and x n . We denote the half-plane weight of this event 
by TLij-.kl-.mni and it is given by (33) with N = 3 and s = 1. This formula contains a double contour integral, and the 
contours must not intersect in order to guarantee a solution of the system of null-state PDEs (appendix A) . According 
to the discussion preceding (2), there are C3 = 5 possible boundary arc connectivities, and in each, z may touch any 
one of the three boundary arcs to give a total of fifteen possible one-pinch-point events. 

Now we associate certain choices of integration contours with particular linear combinations of these configurations. 
In the previously considered cases with one screening charge, we noted that a half-plane pinch-point weight with a 
simple contour connecting z with z by crossing a specified interval (x a ,Xf,) corresponds to a specified pinch-point 
event, and now we investigate to what extent this remains true in our present situation with two screening charges. 
We suppose that 71 connects X\,X2, and z. Then X3,...,xq are connected pairwise by the two remaining boundary 
arcs in one of two possible ways. In both cases, topological considerations show that fusion of the bulk two-leg operator 
^i(z) with its image across the intervals (#3, x 4 ) or (#5, Xq) must give rise to a boundary four-leg operator to leading 
order. Hence, we choose the first contour Ti2 to be a simple curve connecting z and z and crossing the real axis only 
through (2:3, X4). A natural choice for the second contour would be the same as for the first but crossing (x 5 ,xq) 
instead, yet this is not allowed because otherwise this contour would intersect T%2 at z and z. We suppose that the 
second contour is [24,2:5] instead. Now, it is easy to show that 



'(« 2 -«i) 8/k n (ui-xi) 



-4/k 



(x 5 



-4/k, 



du2 du\ = 0, 



(71) 



i=l,2 



where u\ is integrated around a simple loop T surrounding [24, x$\ and the ellipsis stands for the rest of the integrand 
in (33) with N = 3 and s = 1. Then this identity (71) allows us to deform T 12 into a simple curve crossing the real 
axis only through (25, xq), so a bulk- image fusion across (x^^Xq) will produce a boundary four- leg operator there as 
well (figure 11). The formula that follows from these contour choices sum over both possible connectivities of the two 
boundary arcs joining X3,... ,xq pairwise. Their relative coefficients can be found in the usual way. Thus, we have 





n n r\ 



FIG. 11: The decomposition of the integral in (72) into a linear combination of Ili2:34:56 and IIi2:36:45. The contour Ti2 
connecting z with z in the top-left illustration can be deformed into the contour that is a vertical reflection of Ti2 across [xs, xe] 
plus the dashed contour in the top-right illustration. According to (71), integration along the dashed contour gives zero. 
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-tfdb^ o q_ _c£±^ o o 










(columns) may be decomposed into linear combinations of the fifteen one-pinch-point densities 



TABLE I: The M x K k 
(rows). (Relative) coefficients for four of the fifteen decompositions are shown in this table. 



from (33) with N = 3 and s = 1 that 

n | z _ f |K/8+8/K-2 



IIl2:34:56 + nHi2, 



36:45 



ip(-4/K, -4/k) 2 V4 - n 2 



Hixj-xtf^Hlz-x^- 8 ^ I I d Ul du 2 N 



x 5 

r 12 Jx A 



(72) 



The ellipsis stands for the rest of the integrand in (33), and the normalization follows from requiring that we recover 
the two-pinch-point weight IIi2:34 with xq <-> x± upon sending 2:5 — >■ 2:4. Although this density is a natural observable, 
the left side is not a single one-pinch-point configuration. Cyclic permutation of the indices generates only five more 
equations involving just twelve of the fifteen possible weights. One of the missing weights is IIi4 : 23 : 56, and the other 
two missing weights are generated by rotating the hexagon. 

In order to isolate all fifteen weights, we pursue our second strategy of splicing charge-neutral pairs into simpler cor- 
relation functions. We begin with the one-pinch-point weight n 12 (a; i , Xj\ z, z), given by (ipi(xi)^ipi(xj)^i(z)[ 2 ^i(z)). 
We insert into its chiral representation a first charge-neutral collection f*' V-^ 2 (xk)V{^ 2 (xi) V-{u\) du\, chosen so that 
Xk and xi are not separated within the real axis by Xi and Xj, and then we insert a second charge- neutral collection 
fx™ Vi2( x m)V{~ 2 (x n )V-{u2) du 2 , chosen so that x m and x n are not separated within the (one-point compactified) real 
axis by Xi, Xj,Xk, or xi. Wc find fifteen distinct conformal blocks, each with n = m + 1 and I = k + l (i.e., the inserted 
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boundary arcs are not nested) or I = k + 3 (i.e., the inserted boundary arcs are nested) (The case n = 7 or I = 7, 8, 9 
is identified with n. = lorZ = l,2, 3 respectively). Each block has the form 

('0i(^- i )[2]V , i(^j)V'i(a;fc)[o]'0i(2;OV'i(^m)[o]^i(a;n)*i(2 : )[2]*i(2)) = n 2 M(xi, ...,x e ;z, z)Kki;mn{xu ■ ■ ■ , x 6 ; z, z), (73) 

where ICki-mn is the real- valued integral 



fCkl;mn ■= ^(-4//t, -4//t)~ 



du\ du 2 M 



{u x - zfl«-\u x -zfl«- x 

6 

x (u 2 - zf' K -\u 2 - zf /K -\u 2 - Ul f' K l[(ux - Xi )- 4 / K (u 2 - Xl )- 4 ^ 

and where M is given by (33): 

M(xx,...,x 6 ;z,z] 



|: _ 2r /8+8/«-2-Q (x ._ a; . ) 2/«-Q^_ x . 



1-8/k 



(74) 



(75) 



The K-ki-mn can be decomposed into linear combinations of the fifteen crossing weights in the usual way. This 
decomposition is shown for four of the JCki.mn in the top row of table I, and the other eleven are found by cyclically 
permuting the indices. We thus find an invertible system of fifteen equations with the fifteen weights as unknowns. 
The formulas that follow from this inversion are complicated, and we leave their further investigation to the interested 
reader. 



B. Half-plane universal partition functions 



In this section, we construct half-plane universal partition functions from the weights computed in section III A. We 
complete these calculations only for the cases N = 2 and 3. However, the method is clearly generalizable to polygons 
with more sides. In this section, our figures and some of our language suggest that we have conformally mapped the 
upper half-plane onto the interior of a 2.ZV-sided polygon V with the vertices numbered counterclockwise in ascending 
order starting with the bottom-left vertex, with the i-th vertex Wi the image of X4, with xi < Xj whenever i < j, and 
with the bottom side of V sitting flush against the real axis. This language streamlines the discussion of this section, 
and the implementation of this transformation is postponed to the next section IIIC. 

We begin with a background discussion. In the present context, the Kac operator ^i(xj) bears two complementary 
interpretations. First, it is a boundary one-leg operator because it conditions the system so that a single boundary 
arc anchors to the point Xi on the real axis. As described in section II, the 2N boundary one-leg operators in (7) 
generate N non-crossing boundary arcs that cross the upper half-plane and join the Xi pairwise in one of Cjq possible 
connectivities, with Cjv given by (2). If A is an event in which the boundary arcs join in a specified connectivity, then 
we call the (continuum) O(n) model partition function summing exclusively over all samples in A and with boundary 
arcs having fugacity one the "weight" Ii\ of the event A. Because only closed loops should enjoy the full loop fugacity 
n and the boundary arcs do not close into loops (yet), we have endowed them with fugacity one. 

Because the bulk loops have fugacity n but the boundary arcs only have fugacity one, the weight Tl\ is almost, but 
not quite, a physical continuum O(n) partition function. Indeed, we must account for how we physically condition 
the boundary arcs to anchor to the specified boundary points, and this matter will endow the boundary arcs with the 
proper fugacity n. For instance, in the dilute 2 < Q < 4 Potts model, boundary arcs are scaling limits of interfaces 
between clusters with all spins in, say, state A and clusters with all spins in any state but A, so we condition a 
boundary arc to anchor to Xi by changing the be there from wired to spin A to freely assuming all spins except A. 
And in the dense 1 < Q < 4 Potts model, boundary arcs are scaling limits of FK-cluster perimeters, so we condition 
a boundary arc to anchor to Xi by changing the be there from wired to spin A to free. Thus, we find a second 
interpretation of il>x{xi) as a bec operator because it changes the be at Xi from fixed to free or vice versa. (In our 
application, ip\ will also sum over the Q possible spin types for the adjacent fixed segment whenever we are working 
with the Q-state Potts model.) 

In our application, each vertex Wi of the polygon (or in the present setting, its half-plane conformal image Xi) under 
consideration will host a free-to-fixed or fixed-to-free bcc. This sets up an ffbc, of which there are many to consider. 
For example, we can condition a pair of fixed sides to be independently wired, that is, not constrained to exhibit the 
same state, or to be mutually wired, that is, constrained to exhibit the same state. By taking different combinations 
of these options, we generate many different ffbc events. We denote a specified ffbc event by 
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FIG. 12: A boundary loop (red) surrounding an FK boundary cluster. If all sites inside (resp. outside) of the red curve are 
(resp. are not) in spin state A, then the green loop of activated dual bonds is a boundary loop surrounding a spin boundary 
cluster. 



Because its boundary arcs have fugacity one, the weight H\ is independent of our choice of ffbc event ?. To promote 
H\ into a true O(n) partition function that sums exclusively over Af~K and endows the boundary arcs with fugacity n, 
we close the boundary arcs into boundary loops by connecting the vertices of the polygon pairwise via the non-crossing 
exterior arcs mentioned in the discussion preceding (10). (In our present setting, the exterior arcs actually live in 
the lower half-plane.) How this is done is described in [16], and we summarize the details for the rectangle and the 
hexagon. We connect the endpoints of a fixed segment that is independently wired with a single exterior arc, and we 
connect each endpoint of every segment in a collection of mutually wired segments to another endpoint of another 
segment in that collection via an exterior arc. The exterior arcs close the boundary arcs into £ {1, . . . , N} distinct 
boundary loops with each boundary loop contributing a fugacity factor of n (figure 12). Thus, the universal partition 
function summing exclusively over the event A PI <r is n lx *Tl\. By summing over either all boundary arc connectivity 
events A or only those A with samples in the s-pinch-point event A, we find the respective half-plane universal partition 
functions 

T S = ^V^II A , T (Ak) = n ' A " n A- ( 76 ) 

A AnA^0 

There are CV and C s terms in the left and right sum of (76) respectively, with Cjy and C s the N-th and s-th Catalan 
numbers respectively (2). By factoring the pinch-point weight Ha out of the right sum, we obtain (10). 

First, we sum (76) for pinch-point events in the rectangle TZ with the left /right sides wired. When N — 2, there 
are two possible boundary arc connectivity events, enumerated according to figure 13, and there are two possible ffbc 
events, also enumerated according to figure 13 after reflecting the interior arcs outwards into exterior arcs. Thus, <Ti 
is the independent wiring event and q is the mutual wiring event. Both ffbc events are possible in the dense phase 
of the Potts model while only the latter is possible in the dilute phase. If A is a specified one-pinch-point event, then 
the boundary arcs connect in exactly one way, so there is only one term in the sum (76). Using figure 14, we find 
that in the independent wiring ffbc event c^, 

Y(12:34|l) = ™ni2 : 34, T( 3 4 :12 |i) = ^34^2, (77) 




1 2 3 4 5 



FIG. 13: An illustration of the labeling that we will use for the boundary arc connectivities of TZ and Ti. 
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FIG. 14: The four boundary loop configurations for 1Z. Each boundary loop contributes a factor of n. 



T 



(41:23|1) 



n 2 n 



41:23, 



(23:41|1) 



n 2 U 



23:41, 



and we find that in the mutual wiring ffbc event ft, 

Y (12:34|2) = "' 2 ni2:34, 
T(41:23|2) — "-Il4i : 23, 



(34:12|2) 



n 2 n 



34:12, 



T 



(23:41|2) 



112) = nil: 



23:41, 



(78) 

(79) 
(80) 



with the one-pinch-point half-plane weight IlijM an d the loop fugacity n given in (49) and (9) respectively. If A is 
the two-pinch-point event, then for the mutual wiring and independent wiring events, we find 



T(1234|1) — T( 12 34|2) — (" + " 2 )IIi234, 



(81) 



with the two-pinch-point half-plane weight II1234 and the loop fugacity n given in (35) and (9) respectively. 

The half-plane universal partition function T s summing exclusively over the ffbc event <r can be computed from 
(76) by using appropriate generalizations of Cardy's formula for horizontal and vertical crossings in the rectangle [25]. 
The boundary arc connectivity events are labeled in figure 13, so with Tj := T fj and H : 

Ti=nni+n 2 n 2 , T 2 =n 2 IIi+nn 2 

where the weights III and II2 are given by 



n,\ . , we find 



(82) 



r(8/«)r(8/« 

and where G% and G 2 are given by [8] 



1) 



\(x 3 - Xl)(Xi - X 2 ) 



Gi(t?) = G 2 (1 



(l-v)=V 2/K (l-v) 1 - 6/ ^F 1 f-,1--;- V ) 



(84) 



with 2F1 the Gauss hypergeometric function. 

Next, we sum (76) for pinch-point events in the hexagon W with the bottom and top left/right sides wired. When 
N = 3, there are five possible ffbc events. They are enumerated according to figure 13, so ft is the independent wiring 
event, s 2 is the mutual wiring event, and 91 (resp. ci, resp. ft) is the event with the bottom and top-left (resp. bottom 
and top-right, resp. top-left and top-right) sides mutually wired and the remaining fixed side independently wired. We 
call any of the events ?i,?4, and ft a mixed ffbc, and as usual, only the mutual wiring event is possible in the dilute 
phase. If A is a specified one-pinch-point event, then there is only one term in the sum (76). With p the number of 
boundary loops in each sample of A n <r (figure 15), we find 



(ij\kl:mn\s) 



ij:kl:mn j 



pe {1,2,3}, 



(85) 



with the loop fugacity n given by (9) and with the half-plane weight Hij-.ki-.mn given by inverting the system of 
equations that is partly shown in table I. Next, if A is a specified two-pinch-point event, then there arc two terms in 
the sum (76) since the two boundary arcs touching at the two-pinch point can be separated into two nonintersecting 
boundary arcs in two ways. The first and second way has p\ and p 2 boundary loops respectively, with {pi,j>2} = {1,2} 
or {2,3}. Therefore, 



1 {ijkkmnk) = ( nPl + n P2 )Ilijki:: 



{^2} = {1,2} or {2,3}, 



(86) 



with the two-pinch-point half-plane weight Hijki-.mn an d the loop fugacity n given by (70) and (9) respectively. Finally, 
if A is the three-pinch-point event, then there are five terms in the sum (76) since the three boundary arcs touching 
at the three-pinch point can be separated into any of the five possible connectivities shown in figure 13. If <; is either 
the independent wiring ffbc ft or the mutual wiring ffbc ft, then we have 



T (i23456|2) = T (1 23456|3) = ( n + 3" 2 + n 3 )IIi23456, (independent or mutual wiring ffbc) 



(87) 
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FIG. 15: Three possible exterior arc connectivities of H out of five possibilities. Left to right, there are two, two, and three 
boundary loops, giving rise to boundary loop fugacity factors of n 2 ,n 2 , and n 3 respectively. 

and if c is any of the mixed ffbcs <;i, <^4, or <; 5 , then we have 

T (i23456|i) = T (123456|4) = T (1 2 34 56|5) = (2n + 2n 2 + n 3 )ili 23 456, (mixed ffbc) (88) 

with the three-pinch-point half-plane weight II123456 and the loop fugacity n given in (36) and (9) respectively. 

There are many combinations of pinch-point events A and ffbcs events <r for the hexagon. We give explicit formulas 
for some type- A pinch-point densities with the independent wiring event ^3- Except for the three-pinch-point density, 
which is too rare to accurately measure, these results are verified via simulation in section IV. First, the following 
combination sums exclusively over all samples with a specified one-pinch point on a boundary arc connecting vertex 
one with vertex two: 

T(i2:34:56|3) + ™ 2 Y(i2:36:45|3) = « 3 [ni2:34:56 + "1112:36:45], (all fixed sides independently wired). (89) 

The linear combination on the left side is chosen so that the right side equals n 3 times (72). Next, the universal 
partition function for a two-pinch point between two boundary arcs that connect vertices Wq, Wi, W2, and W3 of the 
hexagon is given by (86). We find 

T( 612 3:45|3) = (n + n 2 )Ii6i23:45, (all fixed sides independently wired) (90) 

with n and ILji23:45 given in (9) and (70) respectively. Finally, the three-pinch-point partition function T( 12 3456|3) is 
already given in (87). 

To finish, we give the half-plane universal partition function Tj := T ft summing exclusively over the ffbc event q. 
This is done as follows [19]. If we label the points x\, . . . , xq by a, b, c, d, e, and / in any way with xq — e or /, then 
the ffbc event q with its exterior arcs pairing x\, . . . ,xq into {a, &}, {c, <i}, {e, /} is given by 



Ti(xi, ...,x 6 ) = n 3 ^(~4/ K , -4/«)- 2 Y[(x k - x 3 ) 2 l K ]J(x 6 - x,) 1 ^^ 



j<k 



i=i 



b r d 



3=1 



du\du 2 - (91) 



As usual, Af orders the differences in the integrand so that the integral is real. If k < 4, then we replace the simple 
integration contours with Pochhammer contours that entwine the endpoints and we divide by an extra factor of 
4sin 2 (47r/K). 



C. Transforming the universal partition functions 

In this section, we transform the half-plane universal partition functions to universal partition functions for system 
in the appropriate 2A^-sided polygon. 

We let -Z(A|<r) (resp. Z^) be the half-plane partition function summing exclusively over the pinch-point and ffbc event 
An? (resp. ffbc event where the j-th bec occurs within a small distance Cj from the specified point Xj on the 
real axis in the ffbc event Now, if we send the upper half-plane onto a simply connected domain T> with a smooth 
boundary via a conformal bijection /, then conformal invariance says that 

Zf K ^{wi, . . . ,w 2 n;w) = Z( A |f)(afi, • ■ -,x 2N ; z), Zf(wi, . . .,w 2N ) = Z q (xx, .. .,x 2 n), (92) 
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where ZH (resp. Zf ) is the partition function summing exclusively over the event An<; (resp. <r) in the image system, 
and where Wj — f(xj) € dT> and w — f(x). As the half-plane partition functions £(A|c) an d ^ have the respective 
asymptotic behaviors 

Z (m /Z f ~ C^? 1 . . . £ ^e 2e =T (AK) , V^/ ~ n C^^e? 1 . ..e^T s> (93) 

the image partition functions Zf^) and -Z^ 5 have the respective asymptotic behaviors 

Z^/Z, ^C^iM* ...^W^H' 9 '!^,, ^cf^iM* 1 ...S 2N (w 2N )^Tf , (94) 

with Sj(wj) = €j\df(xj)\, with <5(iu) = e\df(x)\, and with * and respectively given by the usual conformal 
covariance transformation laws 

Tf A|?) := \df(x 1 )\^ . . . \df(x 2N )\- e i\df(z)\- 26 °r (m , Tf := \df(x 1 )\^ . . . \df(x 2N )\- e ^. (95) 

Although it is conformally invariant, the type-A pinch-point density given by the ratio Z'R.JZ'^' is unnatural because 
the radii of the disks containing either the pinch-point or the bees varies with their locations in the closure of T>. For 
this reason, it is natural to replace 6(w) and each 8j(wj) with small numbers S and Sj that are independent of w and 
Wj respectively. The type-A pinch-point density that results is conformally covariant instead of conformally invariant. 

Now we let V be an equiangular 2A^-sided polygon V with vertices at w\, . . . , w 2 n- This is the standard setting for 
this article. To replace S(w) with 5 as prescribed above is valid as long as w is sufficiently far from the vertices. But 
to replace the other 8j{vjj) with Sj is not valid since the derivative of / blows up at each vertex. Instead, we replace 

Z (A|?) ► ^(A|?) g , ^ Q °s C l °1 ■■■°2N° (A|?) ' Z c ' 5 ^ Q C 1 °1 ■■■°2N° L s> l yD J 

where T(a|?) an< ^ ^? are ^ ne respective correlation functions 

Tf A |0 = WW • ■ ■ = . . . i>i(w 2N )) v , (97) 

that use corner one-leg operators in place of boundary one- leg operators at the vertices of V. We define the corner 
one-leg operator by 

i>t{wj) = lim ( 1 ^(wj + £j e ia ), (98) 

£j— >0 y 7T J J 

where <f> is the interior angle of V at the vertex Wj, and where £je la is a complex number (with Ej £ R) such that 
Wj + £je la is on a side of T 3 and very close to Wj. Because V is an equiangular 27V-sided polygon, <j> = (N — l)w/N. 
It is easy to check that if we use corner one-leg operators in place of ordinary boundary one-leg operators, then our 
redefined partition functions (96) are finite. 

First, we transform the N = 2 half-plane universal partition functions into universal partition functions for the 
system in the rectangle 1Z (figure 16). The rectangle has aspect ratio R, vertices w\ = 0,w 2 = R,w% = R + i, and 
ii>4 = i, and interior angles (f> = tt/2. We therefore have 

T(Ak) = (^(OW(if)lif(fl + »)#)$,(M))R = ]im Q (l6eis 2 e 3 e 4 )-^ 

x {M s i)M R - s 2)M R + 'i- s 3)M i + £ ^ s {z,z)) R . (99) 

We use the Schwarz-Christofiel map that conformally sends the upper half-plane onto 1Z and whose continuous 
extension to the real axis sends the points {0, m, 1, oo} counterclockwise to the respective vertices {w%, w 2 , W3, W4}: 

fW = ^T^ A" 1/2 ( m -0" 1/2 (W)- 1/2 dC, r 1 (w) = m S n(wK'\mf. (100) 
2K'(m) J 

Here, K(m) is the complete elliptic integral of the first kind, K'(m) := K(l — m), and m £ (0,1) is the modular 
parameter of the transformation, related one-to-one with R £ (0, 00) through R = K(m) / K 1 (m). This transformation 
(and its continuous extension to the real axis) is invertible, so we define x'j to be the half-plane pre-image of the j-th 
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vertex of 1Z shifted by a small amount £j as shown in (99). These points are given in [26] to leading order, and they 
are respectively near zero, m, one, and infinity. Conformal covariance gives 

4 

r™(m;w,w) = \d.f(z)\- 2e ° hm (16e ie2 e 3 e4)- ei H Idfix^T^iw^w^w^w^w^), (101) 

3=1 

with Tf(A| s ) given in section IIIB. The right side may be computed explicitly, and we find (with z = f (w)) that 

Tf m (m; to, to) = \2mK'sn(wK' \ m)cn(wK' \ m)dn(wK' | m )|[ 16 * 2 -(«-4) 2 ]/8 K 

x [m(l - m)} 6 ^- 1 ^/^ lim x 6/K - 1 T (Ak) (0,m,l,x; z,z). (102) 

x— >oo 

We note that r) — > m as £i, . . . , £4 — > 0. The universal partition function T^? for 7?. is found from (82) in the same 
way after dropping all covariance factors associated with the bulk point w. The result is [16] 



} = n^(mr^ { J£ - m H, F(m) := 2 F, 2 - 12 , 1 - 4 ; 2 - 8 
T^(m) J v y 1 -F(m) J vy \ « « k 



(103) 



Next, we transform the N = 3 half-plane universal partition functions into universal partition functions for the 
system in the hexagon % with vertices W\ — 0,u>2 > 0, 103, . . . , ti)6 in the upper half-plane, and interior angles 
= 271-/3. We have 



T|l k ) = (^ C K)V'l(^2)^(^3)V'l(^4)V'lH5)V'lK)* S (^,w))w 

hm ^ 729^1^4^ ^^ l^ i{w . +E , e U-Ws )Mw ^ . (104) 



We use the Schwarz-Christoffel map that sends the upper half-plane conformally onto % and whose continuous exten- 
sion to the real axis sends the points {0, mi, m2, m 3 , 1, 00} clockwise to the respective vertices {w\ = 0, W2, ■ ■ ■ , w§\ 
of % (figure 16): 

/M = o / C 1/3 (mi - 0" 1/3 (m 2 - C)~ 1/3 (™ 3 - cr 1/3 (i - C)~ 1/3 dC (105) 



3 7o 

This transformation (and its extension to the real axis) is invertible, so we define the point x'j to be the half-plane 
pre- image of the j-th vertex shifted by EjeV~ x > m l 3 . To leading order in £,-, these points are 
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i(l-ma)ea /a , (108) 
>(l-m 3 )e 3/2 , (109) 
)£ 3/2 , (110) 

(111) 



and they are respectively near zero, mi, 771,2, one, and infinity. Conformal covariance gives 



T^| s) (mi,m 2 ,m 3 ;-u;,7«) = |<9/(z)| 26s lim 



£i£2£ 3 £4£s£6 \ 



V 64 / 

6 

xni 9 /(4)r 9lT (Ak)W. w >>>5,»6;»,«') I (112) 

with T(A^) given in section IIIB. The right side may be computed explicitly, and we find (with z = f (w)) that 
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FIG. 16: The transformation of the upper half-plane to the interior of the rectangle and the hexagon and our enumeration of 
the vertices of either polygon. 



Tft| s) (mi, m 2) m 3 ; w, w) = \27z( mi - z)(m 2 - z)(m 3 - z)(l - z )/ 8 |[ 16s M«-4) 2 ]/24 K 

x [mim 2 m 3 (m 2 — TOi)(m 3 — mi)(m3 — m 2 )(l — mi)(l — m 2 )(l — ™3)]' 6 ~ k " 2k 
x lim a; 6 / K_1 T(A|<;)(0, mi, m 2 ,m 3 , l,x; z, z). 



(113) 



We note that (r?,r, cr) — > (mi, m 2 ,m3) as £i, . . . , £g — > 0. The universal partition function for % is found from 
(91) in the same way after dropping all covariance factors associated with the bulk point w. The result is 



T^(mi,m 2 ,m 3 ) = [mim 2 m 3 (m 2 - mi)(m 3 - mi)(m 3 - m 2 )(l - mi)(l - m 2 )(l - m 3 )} (10 k ^ 2k 



x n J /3(-4/K,-4/K)^ 



(«2 - ^i) s/k n^ 4/K (i - ^r 4/K fiK - Ui )^ 4/k 

i=l 3=1 



du2<lui, (114) 



where {a, 6} is the pair of finite endpoints for one of the two exterior arcs and {c, d} is that for the other exterior arc. 
These endpoints are among {0,mi,m 2 ,m 3 , 1}. In the event of the ffbc with all fixed sides independently wired (^ 3 ) 
(91), we have a = 0, b = mi, c = m 2 , and d — m 3 . 



D. Pinch-point densities in polygons 

In this section, we give explicit expressions for some of these densities for the rectangle 7Z (N — 2) and for the 
hexagon H (N = 3). The behavior of the type- A pinch-point density /0(yu<r) m ^ ne P 0i yg° n T 3 * s gi ven by (4) with the 
bulk s-leg exponent Q s given by (12) and with the universal partition functions T^, i and T s computed in section 
III C for the rectangle and the hexagon. 




FIG. 17: Illustration of the type-(12 : 34) pinch-point configuration in the rectangle. The boundary cluster is shaded gray. 
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1. A one-pinch-point density for the rectangle 



The density of type-(12 : 34) one-pinch points in the rectangle 1Z with the left/right sides independently wired 
behaves as (figure 17) 



(1 - to) 6 /*" 1 



2^(2-^,1 

lm[sn(w K' | to) 2 ] 



K'\sn(wK' | m)cn(wK' \ m)&n(wK' | m)| 



\ 1 1 - m) 

k/8-1 



2G 2 + (n 2 - 2)G 4 - nGi - nG 3 
n(n 2 — 4) 



(to, z, z) . (115) 



The Gi are given in (46-46), and n is given in (9). The density p^ 2 -34|2) wrtn tne left/right sides mutually wired is 
found by replacing the argument of the hypergeometric function with 1 — to in (115) and multiplying by n. Figure 
18 shows a contour plot of P(f 2 -34li) ^ or I sm S boundary clusters (k — 16/3). 




0.0 0.5 1.0 1.5 2.0 

FIG. 18: Contour plot of Ising FK cluster (k = 16/3) one-pinch-point density P(i2-34|i)- 



2. A two-pinch-point density for the rectangle 

The density of two-pinch points in the rectangle 1Z with the left /right sides independently wired behaves as (figure 
19) 



P(1234|l)l m > u, ' U; J A -~ °2° 



2 24/ K -2^/6/«- K /8+l( n -l + _ m )]8/K-I 

2^i(2-f,i-^;2-f|i-TO) 

Im [sn(u>-ftT' | to) z 



|sn(wi4T' | m)cn(wi ; C' | m)An{wK' | to) 



k/8+18/k-3 



(116) 




FIG. 19: Illustration of the 



two-pinch-point configuration in the rectangle. Boundary clusters are shaded gray. 
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The density P(^ 2 -34|2) with the left/right sides mutually wired is found by replacing the argument of the hypergeometric 
function with 1 — m in (116). Figure 20 shows a contour plot of PQ234I1) f° r I sm g FrC boundary clusters (k — 16/3). 




FIG. 20: Contour plot of Ising FK cluster (k = 16/3) two-pinch-point density P(i234|i) 



3. A combination of one-pinch-point densities for a hexagon 

Equation (89) is a combination of one-pinch-point densities in the hexagon % with the bottom and top-left /right 
sides independently wired (figure 21). We have 

[P(12:34:56|3) + ^P(12:36:45|3)] (™1 > ™*> ™*\ ™> ®) ,~ CfS^U^ n 2 )" 1 ^* - z\"/***/-* 

x \27z{ mi - z)(m 2 - z)(m 3 - z)(l z)/8^ M Z L^j . (n 7 ) 

i2(mi,rn 2 ,m 3 ) 

Here, is the definite integral 



Ii(m 1 ,m 2 ,m 3 ;z,z) := Af 



i=l,2 

x(« 2 - W i) s/K n ^ 4/K a 

i=l,2 



IK' 



u ? r 4/K TT(m J -uO- 4/K 




+ rT 6 




FIG. 21: Illustration of the type-(12 : 34 : 56) one-pinch-point configuration plus n 2 times the (12 : 36 : 45) one-pinch-point 
configuration in the hexagon. Boundary clusters are shaded gray. 



2G 



with the contour T starting at z, crossing the real axis through either (m 2 , m 3 ) or (1, 00), and ending at z. One can 
show that 1\ is purely imaginary. Also, I 2 is 



/■mi rm 3 

^(mi,ro2,m 3 ) := / / M 

JO J m-2 



n< 4/K (i-^- 4/K n ( 



rrij — Ui 



-4/« 



i=l,2 



.7 = 1,2,3 



(u 2 -ui) 8/K dw 2 dMi. (119) 





-1.0 -0.5 0.0 0.5 1.0 



FIG. 22: Illustration of the type-(6123 : 45) two-pinch-point configuration in the hexagon (left, boundary clusters are shaded 
gray), and a contour plot of the Ising FK cluster (« = 16/3) two-pinch-point density P(l23456|3) (right). 



4- A two-pinch-point density for a hexagon 

The density of two-pinch points in the hexagon T-L with the top and bottom- left/right sides independently wired 
behaves as (figure 22) 

P% 123 . A5 \ 3) (m 1 ,m 2 ,m 3 ;w,w) ~ Cftf 6 /"""/^^ 1 + 1)/3(-4/k, -4/k) 2 \z - z\«l*-*l"-^h{m x , m 2 , mg)" 1 



[mi(m 3 - m 2 )] 1 - 6/K \27z(m 1 - z){m 2 - z)(m 3 - z){\ - z )/ 8 |V3-«/ 24 + 2 A 
[mim 2 m 3 (m 2 — mi){m 3 — mi)(m 3 — m 2 )(l — mi)(l — 7712) (1 — m 3 )] 2 / K 
(2 - n 2 )(H e + H 4 ) + nHi - 2H 2 + nH 3 + n(n 2 - 3)H 5 
(n 2 -4)(ti 2 - 1) 



(120) 



(mi,m 2 ,m 3 ;z, z), 



where J 2 is given in (119) and the Hi are given in (59, 60-65). Figure 22 shows a contour plot of /9?gi23-45|3) f° r I sm g 
FK boundary clusters (k = 16/3). 



5. A three-pinch-point density for a hexagon 

The density of three-pinch points in the hexagon % with the top and bottom- left/right sides independently wired 
behaves as (figure 23) 

pf 6123 . A5l3) (rni, m2 ,rn 3 ;w,w) ~ q C 2 ^/ K -/8 + i (27/8) i6/3 K - K /24-i/ 3( ^ 2 + 3n -i + 1} 



X \ Z {m X - z)(m 2 ~ Z){m 3 - z){\ - z )|4/3-k/24-32/3k 

x |z-zr/ 8 + 32 / K - 4 /,(-4/ K ,-4/ K ) 2 / 2 (m 1 ,m 2 ,m 3 )- 1 , ^ 
where I 2 is given in (119). Figure 23 shows a contour plot of / c, (^ 12 3-45|3) f° r I sm g FK boundary clusters (k = 16/3). 
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IV. SIMULATION RESULTS 



In this section, we present simulation results that verify our predicted pinch-point densities P(i 2 -34|i) anc ^ <°(i234|i) 
(115, 116) for the rectangle and P^ 2 -34-56|3) + n2 <°(i2-36-45|3) an( ^ '°(6i23-45|3) (H^, 120) for the hexagon. We simulated 
critical bond percolation on a 2000x1000 square lattice (aspect ratio R = 2) in a rectangle TZ and critical site 
percolation on a triangular lattice in a regular hexagon Ti inscribed in a 2000 x 2000 rhombus. (Actually, only half of 
the sites in TZ belonged to the physical lattice, and the other sites belonged to the medial lattice. This was for the 
purpose of performing a hull-walk later. This was also true when we sampled Ising FK clusters in H, but it was not 
true when we sampled site percolation clusters in Ti. In the latter case, all of the sites in our simulation belonged to 
the physical lattice.) Using the Swendsen-Wang (SW) algorithm [27], we also sampled critical Ising FK clusters on 
the same lattices in 7Z and Ti respectively. We independently wired the left /right sides of TZ, and we independently 
wired the bottom and top left/right sides of Ti, leaving all other sides free in both situations. Our simulation results 
show good agreement with our predictions. We used Mathematica (version 8) to evaluate all definite integrals that 
appear in our analytic predictions via the "NIntegrate" function, and we used the Schwarz-Christoffel toolbox [28] for 
MATLAB to numerically perform the transformation from the upper half-plane to the interior of Ti. 

To approximate the continuum limit of these critical models, we used very large lattices in our simulations. This 
suppressed the frequency of each s-pinch-point event at every lattice site, so we generated many samples in order to 
compensate this effect. Overall, about sixteen months of computer time (with about a 2 GHz processor) were used 
to sample percolation clusters in TZ and Ti and Ising FK clusters in 71, and over thirty-two months of computer time 
were used to sample Ising FK clusters in Ti, running simulations on many processors simultaneously. 

As usual, we number the vertices of TZ and Ti counterclockwise, starting with the left vertex of the bottom side of 
either polygon as vertex one. In this section, we also shift the hexagon so that its center coincides with the origin. 



A. The rectangle 

To measure the density of percolation pinch points (k — 6), we simulated critical bond percolation on the square 
lattice in TZ. One-pinch points are bonds inside TZ and on the perimeter of a boundary cluster, and two-pinch points are 
red bonds whose activation or deactivation respectively connects or disconnects the left and right boundary clusters. 
Because pinch points occur only on boundary cluster perimeters, our simulations only sampled these perimeters via 
percolation hull- walks on the medial lattice [29-31]. They did not generate entire percolation bond configurations. 

We performed two hull- walks to generate the two boundary cluster perimeters in TZ (figure 24). The first 
(resp. second) walk starts at the medial lattice site a/2 units below vertex one (resp. above vertex three), where 
a is the lattice spacing. One walk ends at the medial lattice site a/2 units below vertex two, and the other ends at 
a/2 units above vertex four. Each step of the walk is located at the midpoint of a bond, and at each step we decide 
to activate or deactivate that bond with critical bond activation probability p c — 1/2. If the bond is activated, then 
the walk turns right. Otherwise it turns left. Each site that is visited by either hull- walk is a one-pinch point, and 
each site that is visited by both hull-walks is a two-pinch point, or red bond. (Actually, one-pinch-point events and 
two-pinch-point events at z € TZ are mutually exclusive according to our definitions. However, our method counts a 



2cS 




FIG. 24: Typical percolation (left) and Ising FK (right) cluster (black, purple, and red bonds) and hull-walk samples (blue and 
green) generated by our simulations in TZ. Two-pinch points are colored red. We note that for percolation (left), our simulation 
only generates bonds that comprise the boundary clusters' perimeters. 



two-pinch point as one-pinch points of both hull walks too. This miscounting is insignificant though because the ratio 
of the two-pinch-point density to any one-pinch-point density vanishes as the system size increases, or equivalently, 
as the lattice spacing a goes to zero, with power a 2 (°2- e i) — a 6 / K .) 

When both of the hull-walks are finished, we note how they connect the vertices of TZ and bin one-pinch-point 
events accordingly. If they connect vertex one with vertex two and vertex three with vertex four (resp. vertex one 
with vertex four and vertex two with vertex three) to create a horizontal (resp. vertical) crossing, then one-pinch-point 
events on the first and second hull- walks contribute to p^ 2 -34|i) anc ^ /°(34-i2|i) ( res P- /°(4i-23|i) an( ^ *°(23-4i|i)) respectively. 
Two-pinch-point events did not need to be sorted this way, although it is interesting to note that the red bonds are 
activated (resp. deactivated) in the event of a horizontal (resp. vertical) crossing. After generating about 2 x 10 s 
samples, we tallied the number of each kind of pinch-point event at each site in an array and divided the total by the 
corresponding array value at the center of TZ to eliminate the non-universal constant and scaling factor appearing in 
(115, 116). With this normalization scheme, our measured densities always equal one at the center of TZ. 

Our percolation pinch-point simulation results are compared with our theory predictions for P(^ 2 -34|i) anc ^ / 9 (i234|i) 
(respectively (115) and (116) with n = 6) in figure 25. In the figure, TZ has been centered at (1,0.5) and rescaled 
to have length two and height one, and we have fixed y to 0.1,0.2,0.3,0.4, and 0.5. The top plot shows the density 
^(i2-34|i) °f P om ts touched by the hull- walk connecting vertices one and two. The left/right peaks of the top curve 
in the figure show that these events are most likely to occur near vertices one and two, and this is expected because 
these vertices are connected by the hull- walk. The bottom plot shows the density /0^ 2 34|i) °f two-pinch points, or red 
bonds. Our results show that these events are most likely to occur near the center of TZ, and this is expected because 
there is more room in the center of 1Z for the hull-walks to collide. Errors averaged over x and the standard deviation 
from this average are shown in table II. 

To measure the density of Ising FK pinch points (k — 16/3), we generated about 4 x 10 8 samples of critical Ising FK 
clusters on the square lattice in TZ via the SW algorithm with critical bond activation probability p c — ^/(l + \/2)- 
In each sample, the left/right sides of TZ were independently wired. That is, all FK bonds of either side were activated 
and necessarily the same color as the other bonds in its boundary cluster, but the bond color of the left boundary 
cluster was allowed to differ from that of the right. The perimeters of the boundary clusters anchored to these 
wired sides were explored via the hull-walk used for percolation to detect the pinch points (without ever activating 
or deactivating FK bonds during the walk) (figure 24). One-pinch-point events and two-pinch-point events were 
detected by the hull-walk in the same way as in bond percolation. However, we note that the "red bond definition" 
of a two-pinch point in percolation does not have a perfect analog with FK clusters. If two boundary FK clusters 
of different colors are separated by a single unactivated bond, we cannot activate this bond and join the clusters, so 
such a bond is not "red" according to the definition of a "spin red bond" in [5]. (However, because the Q colors of 
the Q-state Potts model are distributed uniformly across the FK clusters, the fractal properties of the set of spin red 
bonds and the set of two-pinch-point bonds are the same.) 

Our Ising pinch-point simulation results are compared with our theory predictions for P^ 2 -34| i) an< ^ '°(i234|i) ( rc ~ 
spectively (115) and (116) with k = 16/3) in figure 26. Explanations and notable features of these two plots are 
the same as those for figure 25, except that y is fixed to 0.1, 0.3, 0.4, 0.5, and 0.6. In addition, we note that our 
measurements of the one-pinch-point density consistently exceeded our predictions by a small amount that increased 
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as y decreased. We suspect that this is a finite-size effect that may be reduced by sampling from a larger system. A 
similar effect was observed for the simulations in [16]. 



Avg. error 


y = o.i 


y = 0.2 


y = 0.3 


y = 0.4 


y = 0.5 


y = 0.6 


Perc. s = 


1 


-0.004 


0.023 


-0.003 


0.012 


-0.002 


- 


Perc. s = 


2 


0.010 


0.012 


0.014 


0.014 


0.014 




Ising s = 


1 


-0.080 




-0.062 


-0.048 


-0.036 


-0.346 


Ising s — 


2 


0.008 




0.009 


0.012 


0.014 


0.046 


Std. dev. 




y = 0.1 


2/ = 0.2 


V = 0.3 


y = 0A 


y = 0.5 


y = 0.6 


Perc. s = 


1 


0.002 


0.026 


0.002 


0.025 


0.000 




Perc. s = 


2 


0.003 


0.003 


0.003 


0.005 


0.004 




Ising s = 


1 


0.008 




0.006 


0.013 


0.006 


0.042 


Ising s — 


2 


0.004 




0.006 


0.005 


0.007 


0.034 



TABLE II: The error (theory minus simulation) averaged over x, and the standard deviation of the error from that average, of 
the data displayed in figures 25 and 26. 
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FIG. 25: Density of percolation [k = 6) one-pinch points P(i2-34|i) (top) and two-pinch points P(i234|i) (bottom) versus i in a 
2000x1000 rectangle rescaled to 2 x 1 for various values of y. Both densities are normalized to equal one at the center (1,0.5) 
of TZ. 
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FIG. 26: Density of Ising FK = 16/3) one-pinch points /9(i2-34|i) (top) and two-pinch points P(i234|i) (bottom) versus i in a 
2000x1000 rectangle rescaled to 2 x 1 for various values of y. Both densities are normalized to equal one at the center (1,0.5) 
of TZ. 



B. The hexagon 

To measure the density of percolation pinch points, we generated about 6.56 x 10 s samples of critical site percolation 
on the triangular lattice (site activation probability p c — 1/2) in H. Again, because pinch-point events occur on the 
perimeters of the three boundary clusters, our simulations sampled only these perimeters via three distinct site- 
percolation hull-walks on the triangular lattice (left picture in figure 27). The first, second, and third hull-walks 
respectively started at vertices one, three, and five and in the direction pointing into the adjacent free side of H. 
Each hull-walk necessarily ended at an even vertex, and no two hull-walks could end at the same vertex. The finished 
hull-walk actually consists of two juxtaposed paths of neighboring activated and deactivated sites. We will call the 
former (resp. latter) path the inner (resp. outer) boundary arc of the boundary cluster (figure 27). A path, called a 
smart kinetic walk, passes between these juxtaposed paths with each step on the dual (honeycomb) lattice [32]. 

Pinch-point events were counted in almost the same way as in the simulations for the rectangle. Each point on the 
inner boundary arc is a one-pinch point, activated (resp. deactivated) sites shared between two inner or outer boundary 
arcs are two-pinch points, and activated (resp. deactivated) sites shared between three inner or outer boundary arcs are 
three-pinch points. Again, although one-pinch-point, two-pinch-point, and three-pinch-point events are technically 
mutually exclusive, we counted three-pinch-point events as two-pinch-point events and two-pinch-point events as 
one-pinch-point events in our simulations, knowing that this over-counting has a negligible effect on our results. 
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FIG. 27: Typical percolation (left) and Ising FK (right) cluster samples H generated by our simulations. In the left illustration, 
black, red, blue, and green (resp. gray, yellow, pink, light blue, resp. white) sites are activated (resp. deactivated, resp. undecided), 
black paths are smart kinetic walks, red, blue, and green (resp. gray) sites form the inner (resp. outer) boundary arc of a boundary 
cluster, and yellow, pink, and light-blue sites are two-pinch points. In the right illustration, red sites and centers of red bonds 
are two-pinch points. 



Our percolation pinch-point simulation results are compared with our theory predictions for P(^ 2 -34-56|3) ~'~' (i2-36-45|3) 
and P(6i23-45|3) (respectively (117) and (120) with k = 6) in figure 28. Three-pinch-point events were so rare that 
we could not generate enough samples to verify our three-pinch-point density prediction (121). In the figure, H has 
been centered at the origin, rescaled to have side-length one, and oriented as in figure 27, and we have fixed y to 
—0.69, —0.52, —0.31, and —0.03. The top plot shows the density P^ 2 -34-56|3) ' 9 (i2-36-45|3) 01 one -pinch points touched 
by the hull-walk connecting vertices one and two. (This density includes all such events without regard to how the 
other two walks connect vertices three through six.) As expected, this density is greatest near vertices one and two 
which are connected by the hull-walk. The bottom plot shows the density P(^ 12 3-45|3) °f two-pinch points, or sites 
touched by a hull-walk connecting vertices one and two and another hull walk connecting vertices three and six. This 
density grows at first as we move from the bottom side of H towards its center, as expected, and then it diminishes 
as the center is approached. This diminishing is also expected since P?gi23'45|3) * s greatest below the center of the 
rectangle conformally equivalent to % with vertices one, two, three, and six, and the image of this point in % is below 
the origin and on the y-axis. 

To measure the density of Ising FK pinch points, we generated 5.68 x 10 8 samples of critical Ising FK clusters on 
the triangular lattice in H via the SW algorithm. The critical probability of FK bond activation on the triangular 
lattice is p c = (y/3 — l)/v / 3 f=s 0.42265 [33]. Because the triangular lattice is not self-dual, a hull-walk similar to that 
of 72 that steps on both regular and dual sites equally often is not so easy to implement. For this reason, we chose 
our hull- walk to step along the bonds in % forming the perimeter of a boundary FK cluster (i.e., the inner boundary 
arcs). We also tracked the activated dual bonds that surrounded the boundary clusters (i.e., the outer boundary 
arcs). The s-pinch-point events were identified with intersections between s distinct inner or outer boundary arcs in 
the same way as with percolation (figure 27). 

Our Ising FK cluster pinch-point simulation results are compared with our theory predictions for Pn2-34-56|3) + 
2^(12 -36 -45 1 3) an d P(6i23-45|3) (respectively (117) and (120) with k — 16/3) in figure 29. Explanations and noticeable 
features of these two plots are the same as those for figure 28. The top plot shows the density Pn2-34-56|3) ~'~^' (i2-36-45|3) 
of one-pinch points touched by the hull-walk connecting vertices one and two. The factor of two arises from the factor 
of n 2 in (117) and because n — v2 when n — 16/3. This factor may seem unnatural, but omitting it considerably 
worsens the agreement of prediction and simulation. Thus, this agreement gives a non-trivial verification of our 
prediction for the relative coefficients between the two densities in (117). We note that this one-pinch-point density 
deviates from our prediction by much less than that in the rectangle. This seems peculiar since the system size in 
the rectangle and in the hexagon are about the same. One might suspect that the larger coordination number of the 
triangular lattice over that of the rectangular lattice causes this effect, but this suspicion is contradicted by the large 
errors in the hexagon two-pinch-point case. We therefore do not propose a possible explanation for this observation. 
The deviations of our two-pinch-point simulation results from our prediction (120) with k = 16/3 are larger than 
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^(12:34:56|3) + P(12:36:45 13) 
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FIG. 28: Density of percolation (k — 6) one-pinch points P(i2-34-56|3) + P(i2-36 45|3) (top) an d two-pinch points P(6i23-45|3) 
(bottom) versus i in a regular hexagon inscribed in 2000x2000 rhombus and adjusted to have side-length one and center at 
the origin. Both densities are normalized to one at the center of T-L. 



those for percolation. This may be understood as follows. The two-pinch-point density scales as a 2 as the lattice 
spacing a goes to zero [34]. For percolation (k = 6), the power is 5/4, but for Ising FK clusters (n — 16/3), the power 
increases to 35/24. Thus on a very large lattice, the frequency of Ising FK two-pinch-point events is suppressed even 
more than that of the percolation two-pinch-point events. Indeed, this is what have observed. We expect that more 
samples would lessen the deviation. However, such simulations would require very considerable computer resources. 



V. SUMMARY 



We summarize the main results of this article. The half-plane weight Ha of the type-A s-pinch-point configuration 
is given by (49, 35) for the cases N = 2 and s = 1 and 2 and is given by (72, 70, 36) for the cases N = 3 and s = 1, 2, 
and 3 also for various pinch-point events A. After specifying a particular ffbc event we construct the universal 
partition function T(a| ? ) that sums exclusively over the event Afl;, and we transform it into a universal partition 
function with the rectangle 1Z or the hexagon H for its domain. The results are (102) and (113) respectively. Then 
to within non-universal factors, the type-A pinch-point density P(a|<j) is found by dividing the transformed universal 
partition function by the universal partition function summing exclusively over the ffbc event <j. This prediction 
agrees well with measurements by simulation that sampled one-pinch-point and two-pinch-point events between critical 
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^(12:34:56|3) + z ^{12:36:45 1 3) 



y = -0.69 



y = -0.52 




2/ = -0.31 
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FIG. 29: Density of Ising FK (k = 16/3) one-pinch points P(i2-34-56|3) + 2p(i2-36'45|3) (*°P) an d two-pinch points P(6i23-45|3) 
(bottom) versus i in a regular hexagon inscribed in 2000x2000 rhombus and adjusted to have side-length one and center at 
the origin. Both densities are normalized to one at the center of T-L. 



percolation and Ising FK boundary clusters in rectangles and hexagons with every other side independently wired. 



Avg. error 


y = -0.69 


y = -0.52 


y = -0.31 


y 


= -0.03 


Perc. s = 1 


-0.005 


-0.004 


-0.003 




-0.002 


Perc. s — 2 


-0.017 


-0.024 


-0.024 




-0.013 


Ising s = 1 


0.006 


0.001 


-0.000 




-0.000 


Ising s = 2 


0.140 


-0.002 


0.005 




0.068 


Std. dev. 


y = -0.69 


y = -0.52 


y =-0.31 


y 


= -0.03 


Perc. s — 1 


0.004 


0.004 


0.003 




0.002 


Perc. s — 2 


0.011 


0.015 


0.013 




0.006 


Ising s = 1 


0.006 


0.006 


0.004 




0.002 


Ising s = 2 


0.052 


0.114 


0.107 




0.045 



TABLE III: The error (theory minus simulation) averaged over x, and the standard deviation of the error from that average, 
of the data displayed in figures 28 and 29. 
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Appendix A: A solution space of the system (13-16) 



We demonstrate that the Coulomb gas solution (33) solves the null-state PDEs (13) and the Ward identities (14-16). 
This calculation is modeled after a similar calculation presented in [35]. Throughout this section, we use the dense 
phase (k > 4) conventions for our Kac charge and screening charge notations (29), we modify the notation of (33) 
by letting X2N+1 = z,x 2 n+ 2 = z, and X2N+2+j '■= Uj, and we assume that s < N. When s = N, (33) reduces to 
(31), and because this formula is purely algebraic, it is easy to directly verify that it solves the system. As usual, 

X\ < X 2 < . . . < X2N-1 < %2N- 

To begin, we consider the function 



2N+2+M 



$(xi,...,X 2 N+2+m) = JJ {Xj-Xif 



(Al) 



%<3 

Our strategy is to choose the a, (i.e., the charges of the chiral operators) such that for 1 < i < 2N, 

2N 

i^ + E 

j^i - - 

2N+2+M 



d, 



(Xj Xi)' 



0: 



2N+1 



X2N+1 - X{ 
92N+2 

X2N+2 - Xi (x 2 N+l - Xi) 2 {x 2 N+2 ~ X l ) 2 



E d ^ 

k=2N+3 



(A2) 



where we recognize the differential operator of the z-th null-state PDE (13) on the left side of (A2), and then integrate 
X2N+3, ■ ■ ■ , x 2 n+2+m on both sides of (A2) around nonintersecting closed contours Ti, . . . ,Tm such as Pochhammer 
contours entwining pairs of the points x\, . . . , 2:2^+2 (figure G). Integration on the right side gives zero. On the left 
side, the integrand is a smooth function of X\,... , x 2 n+M+2 because the contours do not intersect. Thus, we may 
commute each differentiation with each integration to find that the M-fold integral § $ solves the system (13). 
With some algebra, we find that for any choice of {hj}, {otj}, M G Z + , and 1 < i < 2N, 



2N+2+M 



E 



o,, 



[Xj X{)' 



2N+2+M 



E 



aja k (na1 - 1) 



2AT+2+M 



E 



ctiajlnaiaj — k/2 + 2) — hj 

(Xj — Xi) 2 



$. (A3) 



If we set hj = 6 X for 1 <j < 2N, hj = O s for j = 2N + 1, 2N + 2, and hj = 1 for j > 2N + 2 (the conformal weight 
of the V± chiral operators that generate screening operators), then for 1 < i < 2N, we can write (A3) as 



2N 

E 

2N+2+M 

E 

3^k 



[Xj Xi)' 



()■ 



2N+1 



2N+2 



e. 



X2N+1 - X l X 2N + 2 - X t (x 2N +\ - Xi) 2 (x 2N+2 



OljOlk^KOt 2 



2N+2+M 



(Xj X^i^Xfc Xi) 



E 

3=£i 



a,ittj{KUiOij — k/2 + 2) — hj 



{Xj Xi)' 



2N+2+M 
k=2N+3 



X k Xi 



Next, we choose a particular i E {1, . . . , 2N}. If we set 



= a 12 = \ j\fn 1 otj — cto ± yi 



a 



'j • 



(A4) 



(A5) 



35 



then the term in brackets on the right side of (A4) vanishes, giving the desired form (A2). In order to satisfy (A5) 
for all 1 < i < 27V, we need 

a j = a^ 2 , l<j< 27V, (A6) 

aj=a^ s , 27V + 1 < j < 27V + 2, (A7) 

aj = a±, 27V + 3 < j < 2N + 2 + M. (A8) 

Explicit formulas for (A6-A8) are given by (29). If we make this choice, then <f $ solves the 27V null-state equations. 

To finish, we choose M and the signs of the square roots in (A7, A8) so that F := § $ solves the conformal Ward 
identities (14-16) too. This happens if the sum of the powers in $ involving the screening charge otj equals negative 
two for each j > 27V + 3. We now motivate this claim. The Ward identities dictate that F must have the conformally 
covariant form (17), which is equivalent to demanding that the function 

2JV-1 

G(xt, ...,x 2N ;x 2 n+i,X2N+2) ■= (x 2 N+2 - x2N+i) 2eB J| {Xj+1 ~ Xj) 26l F(xi, , . . , X2N+2) (A9) 

j — 1 , odd 

2N-1 

= (.x 2W+2 -x 2JV+1 ) [i6s2 - (k - 4)2i/8k n (z, + i-z,) 6/k - 1 (Aio) 

j — 1, odd 

X <j) . . . <j) $(x 1 , . . . ,X 2 N+2+M)dx 2 N+3 ■ ■ ■ dx 2 N+2+M (AH) 

only depends on a set of 27V — 1 independent cross-ratios that can be formed from x±,.. . ,X2N+2- in particular, we 
choose 

V 3 =f(x 3 ), j = 2,..., 27V - 2, 27V + 1,27V + 2, where f(x) := & ~ x i)( x ™ ~ x ™-i) (M2) 

(X2N-1 ~ Xi){X 2N - X) 

and we note that f(x\) — < 772 < 1/3 < . . . < r)2N-2 < /few-i) = 1 < J{x2n) = 00. Then this condition amounts 
to requiring that G satisfy 

G(xi,x 2 ,x 3 , . . . , x 2N - 2 ,X2N-i,X2N] x 2 n+i,x 2 n+2) = G(0, rj 2 ,r] 3 , . . . , r\ 2 N-2i 1, oo; ri2N+i,mN+2)- (A13) 

Because the right side must ultimately be finite (and this requirement fixes M once the signs in (A7, A8) are chosen), 
we momentarily ignore all of its infinite factors. Thus, the j-th integral on the left side of (A13) has the form 

/2AT+2 2W+2+M 
H(x k - Xj f* n (x^x^dx, (aw) 

k=l l=2N+3 

while the j-th integral on the right side of (A13) has the form (using rjj instead of Xj for our symbol of integration) 

„ 2N+2 2N+2+M 

/^(i-^y— J] fa-ntf* J] (m-Vjf'dvj. (A15) 

J k=2 l=2N+3 

k^2N-l,2N l^kj 

(We note that the integrand of (A14) contains an extra factor that was ignored in (A15) when x 2 n was set to infinity.) 
The simplest way to achieve the equality in (A13) is for these two integrals to be the same up to algebraic prefactors. 
Upon introducing the change of variables 

% = /(^)> 3 + 1,27V -1,27V, (A16) 
with / defined in (A12), the integral (A15) transforms into 

V( Xl ,...,x 2N+ 2) II (-^-T-) II f^- 52 -) I 3 (A17) 

J \X2N-XjJ l=2 ^ +3 \X2N-Xj) (x 2N -X 3 y 

k^2N l^j 



3G 



where V(x\, . . . ,X2N+2) is an algebraic prefactor whose explicit formula is presently unimportant. The main point 
for now is that in order to match the integral in (A17) with that in (A14), we must have 

@2N = —Pi — ... — ~ (3j + i — ... — P2N-I — 02N+1 — ... — P2N+M+2 — 2. (A18) 

Thus the sum of the powers in the integrand of the j-th integral (A14) appearing in § $ (when that integral is 
considered to be evaluated first) should equal negative two. This must be true for all j G {2N + 3, . . . , 2N + M}. 
We note that this requirement is satisfied if all of the aj in (AI) sum to 2ao, with 2ao = <x+ + a_ the background 
charge; this is the Coulomb gas neutrality condition [20]: 

- 2 = ^ A = 2oy X! ai = 2a i X! a% ~ a i) = 2a± ai ~ a± ) < * => X! a% = 2a °- ( Al9 ) 

i^j i^j \ i / \ i / i 

(In (A19), we have used the identities a + + a_ = 2a and = —1.) Thus, we should choose M and the signs of 

the square roots in (A7, A8) so that the total charge J2j a j equals 2o.q and § $ therefore solves the Ward identities. 
We suppose that m + (resp. m_) of the screening charges equal a + (resp. aJ) with m + + m_ = M necessarily. Then 
the total charge is 

2ao — aj — 2Na.i 2 + a^ s + a^ s + m+a+ + m-a.- (A20) 
3 

= -Na- + a + + (1 ± s)a_/2 + (1 ± s)a_/2 + m + a + + m-Ot- (A21) 
(m+ + l)a+ + (-N + s + m- + l)a_ ++ 

(m + + l)a + + (-N - s + m_ + l)a_ — . (A22) 
(m+ + l)a+ + (-N + m_ + l)a_ H — 



= < 



The H — h (resp. ) case corresponds to using the + (resp.—) sign for both ao jS , and the H — case corresponds to 

using the + sign for one of the qo,s and the — sign for the other. In order for these expressions to equal 2oq, we need 
the coefficients of a + and a_ in (A22) to equal one. Therefore, to + = and m_ = M, so all M screening charges 
use the — sign. The number of screening charges is thus given by 



M = { N + s . (A23) 




The ++ case is used exclusively in this article since it apparently contains all of the pinch-point densities when 
N = 1,2, and 3. We anticipate that the ++ case contains all pinch-point densities for all N 6 Z + . Moreover, we see 
that these solutions span a subspace of a larger solution space containing all cases in (A23) (among possibly more 
solutions). 

To explicitly show that § $ satisfies the conformal Ward identities (14-16), we can check that relation (A13) is 
satisfied by introducing the substitution rjj = f(xj), with j = 2, . . . , 2A^ — 2, 2A^ + 1, . . . , 2A^ + 2 + M and with / 
defined in (A12), on the right side and performing some straightforward algebra. 
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